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ABSTRACT 


A  study  is  made  of  thin  inextensible  elastic  rods  undergoing 
large  deflections  when  subject  to  concentrated  and  distributed  loads. 

A  general  numerical  method  is  developed  for  solution  of  the 
problem  of  determining  the  deformed  shape  of  a  rod  when  the  loading 
and  boundary  conditions  are  specified.  Some  specific  boundary  value 
problems  are  solved  including  the  problem  of  a  tapered  cantilever 
bent  under  its  own  weight.  Experimental  verification  of  the  numerical 
method  is  given  for  the  problem  of  a  uniform  rod  pinned  at  both  ends 
and  hanging  under  its  own  weight. 

The  inverse  problem  of  determining  the  loading  of  a  rod 
needed  to  deform  it  to  a  specified  shape  is  also  considered.  Some 
special  analytic  solutions  are  presented  and  a  numerical  integration 
method  is  used  for  more  general  shapes. 
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CHAPTER  I 

INTRODUCTION 

In  recent  years  considerable  attention  has  been  focused  on 
the  production  of  oil  and  gas  reservoirs  through  offshore  pipelines. 
The  laying  of  these  pipelines  in  deep  water,  particularly  in  the 
Arctic,  can  be  extremely  difficult.  One  method  of  installing  these 
pipelines  is  to  lay  it  continuously  from  a  barge  or  through  a  trench 
in  the  ice.  To  avoid  buckling  or  overstressing  the  pipe  during  this 
procedure  the  radius  of  curvature  of  the  pipe  must  be  kept  above  some 
minimum  allowable  value.  This  can  be  done  by  supporting  the  pipe  in 
a  bouyant  stinger  or  flexible  truss  as  shown  in  Figure  1.  In  design- 
ing  this  supporting  structure  two  types  of  problems  must  be  solved. 

In  the  upper  supported  region  the  designer  can  choose  an  optimum  shape 
for  the  pipe.  The  supporting  loads  required  to  give  the  pipe  this 
shape  must  then  be  found.  In  the  lower  unsupported  region  the  pipe 
deforms  due  to  its  weight.  The  loading  of  the  pipe  in  this  region 
is  then  known  and  to  check  the  radius  of  curvature  the  deformed  shape 
of  the  pipe  must  be  determined. 

These  problems  can  be  solved  by  treating  the  pipe  as  a  thin 
inextensible  rod  loaded  by  both  distributed  and  concentrated  loads. 

As  in  the  usual  Bernoulli-Euler  beam  theory  the  bending  moment  at  any 
point  along  the  rod  is  proportional  to  the  curvature  at  that  point. 

In  the  small  deflection  theory  this  relationship  is  linearized  by 
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SUPPORTING  VEHICLES 


(b) 


FIGURE  1:  Offshore  Pipeline  Installation 
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neglecting  the  square  and  higher  powers  of  the  slope  of  the  rod  as 
small  compared  to  unity.  For  this  problem  the  square  of  the  slope 
of  the  pipeline  cannot  be  neglected  and  the  nonlinear  relationship 
must  be  used. 

The  general  problem  of  nonlinear  bending  of  thin  rods  requires 
three  basic  components.  These  are  the  free  or  unloaded  shape  of  the 
rod,  the  loading  of  the  rod,  and  the  deformed  shape  of  the  rod.  The 
usual  problem  considered  is  that  in  which  the  free  shape  and  loading 
of  the  rod  are  given.  The  deformed  shape  of  the  rod  must  then  be 
determined.  The  inverse  problem  is  that  of  finding  the  free  shape 
of  the  rod  when  the  loading  and  the  deformed  shape  are  specified.  A 
further  problem  is  that  in  which  the  free  and  deformed  shapes  of  the 
rod  are  given  and  the  required  loading  is  to  be  found.  The  two 
pipeline  problems  mentioned  previously  are  special  cases  of  the 
first  and  third  problems  given  above  where  the  free  shape  of  the  rod 
is  a  straight  line.  The  second  problem  of  determining  the  free  shape 
of  the  rod  is  not  applicable  to  the  study  of  initially  straight  rods. 

In  Chapter  II  the  governing  differential  equations  for  the 
three  general  problems  given  above  are  derived.  These  equations  are 
valid  for  a  rod  of  varying  cross  section  and  subject  to  nonuniform 
distributed  loading.  Chapter  III  considers  the  problem  of  determining 
the  deformed  shape  of  an  initially  straight  rod  loaded  by  distributed 
and  concentrated  loads.  A  general  numerical  method  is  presented 
which  can  be  easily  applied  to  the  pipeline  problem  and  to  problems 
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of  rods  with  varying  cross  section  and  nonuniform  distributed  loadings. 
With  this  numerical  method  only  a  small  programable  desk  computer  with 
limited  memory  is  needed  to  give  accurate  solutions.  The  numerical 
method  is  compared  to  two  of  the  few  analytical  solutions  known  for 
problems  of  this  type.  Experimental  verification  of  the  numerical 
method  is  presented  in  Chapter  V  for  the  case  of  a  thin  rod  deformed 
under  its  own  weight  when  pinned  at  both  ends.  The  problem  of 
finding  the  distributed  loading  necessary  to  deform  an  initially 
straight  rod  to  specified  shape  is  treated  in  Chapter  IV.  Some 
special  analytic  solutions  are  given  and  a  numerical  integration 
procedure  is  used  to  solve  more  general  problems. 


5. 


CHAPTER  II 

BASIC  CONSIDERATIONS 

2. 1  Classification  of  Previously  Solved  Problems 

The  problem  of  nonlinear  bending  of  thin  rods  has  received 
considerable  attention  in  the  past.  The  usual  problem  solved  is  that 
of  determining  the  deformed  shape  of  the  rod  when  the  free  shape  and 
loading  are  specified.  The  elastica problem  of  an  initially  straight 
rod  bent  in  a  plane  curve  by  couples  and  forces  at  its  end  points  was 
first  solved  by  Euler  [25].  The  problem  of  an  initially  straight 
horizontal  cantilever  with  a  vertical  concentrated  load  at  its  free 
end  was  treated  approximately  by  Boyd  [4]  and  Gross  and  Lehr  [12]. 
Barton  [1]  and  Bisshopp  and  Drucker  [3]  later  solved  this  problem 
in  terms  of  elliptic  integrals.  Analytical  solutions  were  obtained 
by  Conway  [7]  and  Sato  [20]  for  the  case  of  an  initially  curved 
cantilever.  Saelman  [19]  considered  the  problem  of  horizontal  and 
vertical  loads  acting  simultaneously  at  the  free  end  of  an  initially 
straight  cantilever.  The  same  problem  has  been  solved  by  Mitchell  [17] 
for  an  initially  curved  rod.  Frisch-Fay  [9]  considered  the  case  of  a 
cantilever  loaded  by  several  concentrated  vertical  loads.  The 
problem  of  a  simply  supported  rod  with  a  central  concentrated  load 
has  been  solved  by  Conway  [6]  and  Gospodnetic  [11]. 

For  problems  involving  uniformly  distributed  loads  very  few 
analytical  solutions  have  been  obtained.  Hummel  and  Morton  [14] 
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presented  an  approximate  solution  for  the  case  of  a  horizontal 
cantilever  subject  to  a  vertical  uniformly  distributed  load.  This 
same  case  was  treated  by  Biclcley  [2]  and  later  by  Rohde  [IS]  using 
an  approximate  series  solution.  Seames  and  Conway  [22]  produced 
a  tabular  method  using  tangential  circular  arcs  which  could  be  used 
to  approximately  solve  this  problem  and  others,  including  problems 
with  initially  curved  rods.  Iyengar  and  Rao  [15]  used  a  series 
solution  to  solve  the  problem  of  a  simply  supported  rod  loaded 
simultaneously  by  a  uniformly  distributed  load  and  a  central  concentrated 
load.  The  case  of  a  cantilever  subjected  to  a  uniform  normal  load 
has  been  solved  analytically  by  Mitchell  [17].  Lippmann,  Mahrenholtz 
and  Johnson  [16]  have  used  an  analog  computer  to  solve  various 
distributed  load  problems,  including  that  of  a  cantilever  with  a 
varying  distributed  load  and  flexural  rigidity.  The  monograph  by 
Frisch-Fay  [9]  includes  mainly  solutions  of  this  first  type  of  problem. 
Reference  is  made  to  many  of  the  analytical  and  numerical  solutions 
mentioned  above. 

More  recently  various  numerical  methods  have  been  used  to 
solve  these  problems,  some  giving  more  accurate  results  for  problems 
previously  considered.  A  finite  difference  method  has  been  used  by 
Wang,  Lee  and  Zienkiewicz  [29].  Wang  [28]  has  proposed  a  numerical 
method  which  can  be  used  in  problems  with  several  concentrated  loads. 

His  proposed  method  for  rods  with  distributed  loads  is  in  error  as 
pointed  out  by  Holden  [13].  A  finite  element  procedure  has  been 
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developed  by  Tada  and  Lee  [24]  for  nonlinear  bending  problems.  Schmidt 
and  Da  Deppo  [21]  have  used  an  extension  of  Rhodes  [18]  series  solution 
to  provide  a  more  accurate  solution  for  the  deformed  shape  of  a  canti¬ 
lever  with  a  uniformly  distributed  load  and  to  provide  a  solution  for 
the  post  buckling  behavior  of  a  heavy  column.  The  nonlinear  bending 
of  a  cantilever  of  variable  cross  section  has  been  considered  by 
Verma  and  Krishna  Murty  [27].  Recently  a  fourth  order  Runge-Kutta 
method  has  been  used  by  Holden  [13]  to  give  accurate  solutions  for 
problems  with  uniformly  distributed  loads.  Yang  [30]  has  used  a 
matrix  displacement  method  to  find  the  deformed  shapes  of  beams  and 
frames  with  concentrated  and  distributed  loads. 

The  problem  of  finding  the  free  shape  of  a  rod  necessary  to 
give  a  specified  deformed  shape  when  a  known  loading  is  applied  has 
received  some  attention  in  the  past.  Truesdell  [26]  found  the  free 
shape  of  a  cantilever  which  becomes  straight  when  deformed  by  a  uniformly 
distributed  load.  The  general  solution  for  the  free  shape  of  a 
cantilever  with  a  specified  deformed  shape  and  loading  was  given  by 
Mitchell  [17]. 

The  third  type  of  problem,  that  in  which  the  free  and  deformed 
shapes  are  specified  and  the  required  loading  is  to  be  found,  has 
received  little  attention  in  the  past.  This  type  of  problem  occurs  in 
the  bending  of  bandsaw  blades  [10]  and  in  the  bending  of  a  pipe  to  a 
specified  shape  [8]. 
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2 . 2  Basic  Equations 

As  in  all  of  the  investigations  referenced  in  section  2.1 
this  work  considers  the  deformation  of  thin  inextensible  rods  loaded 
in  one  plane  by  static  concentrated  and  distributed  loads.  It  is 
assumed  that  the  loading  and  geometry  of  the  rod  are  such  that  the 
rod  deforms  only  in  the  plane  of  loading  due  to  bending  alone.  The 
rod  is  assumed  to  be  thin  so  that  the  longitudinal  strains  in  the 
rod  remain  small.  The  relative  rotations  in  the  rod  can  still  be 
large  and  this  leads  to  large  deflections  of  the  rod.  According  to 
the  Bernoulli- Euler  law  the  bending  moment  M  at  any  point  in  the 
rod  is  proportional  to  the  change  in  curvature  of  the  rod.  Thus 


M  =  El 


"d0 

ds 


ds 


(2.1) 


where  El  is  the  flexural  rigidity  of  the  cross  section  of  the  rod  and 
<j)  (s)  and  0  (s )  define  respectively  the  free  and  deformed  configurations 
of  the  rod  as  shown  in  Figure  2. 

The  equilibrium  equations  of  the  rod  are  given  by  Frisch- 
Fay  I9J  and  can  be  obtained  from  consideration  of  a  small  section  of 
the  deformed  rod  as  shown  in  Figure  3. 

T,  V  and  M  are  the  tension,  shear  force  and  bending  moment 
respectively.  The  distributed  load  per  unit  length  can  be  resolved 
into  components  w^(s)  and  wy(s)  acting  parallel  to  the  x  and  y  axis 
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FIGURE  2: 
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FIGURE 
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Free  and  Deformed  Rod  Configurations 


T  +  dT 


FIGURE  4:  Section  of  Rod  with  Normal  and  Tangential  Loads 
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or  into  components  wn(s)  and  wt(s)  acting  normal  and  tangential  to 
the  rod  as  shown  in  Figure  4.  The  static  equilibrium  of  forces 
resolved  in  the  x  and  y  directions  requires 

+t£Fy  =  0  =  -  Wyds  +  VcosG  -  (V  +  dV)  cos (9  +  d9) 

-  TsinO  +  (T  +  dT)  sin(0  +  d9)  =  0 


After  expanding  and  neglecting  second  order  terms  this  equation  reduces 
to 


d_ 

ds 


(TsinQ) 


-  (VcosG)  -  w  =  0 
ds  y 


(2.2) 


Similarly 


±£F  =  0  =  w  ds  -  Tcos6  +  (T  +  dT)  cos(0  +  dG) 

X  X 

+  (V  +  dV)  sin (Q  +  dG)  -  Vsin6  =  0 


which  reduces  to 

4"  (TcosG )  +  (VsinQ )  +  w  =  0  (2.3) 

ds  ds  x 

In  some  applications  it  is  more  suitable  to  consider  these  equations 
resolved  in  the  normal  and  tangential  directions.  This  can  be  done 
directly  by  considering  the  static  equilibrium  of  forces  resolved  in 
these  directions  or  by  first  performing  the  indicated  differentiations 
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in  equations  (2.2)  and  (2.3)  and  then  multiplying  equation  (2.2)  by 
cos0  and  equation  (2.3)  by  -  sin0.  Adding  the  resulting  equations 
gives 


d0  dV 

T  ~r~  —  — —  =  w  sin0  +  w  cos0  =  w 
ds  ds  x  y  n 


(2.4) 


Similarly  multiplying  (2.2)  by  sin0  and  (2.3)  by  cos0  and  adding  the 
two  equations  gives 


—  +  V  — 
ds  ds 


-  w  cos0  +  w  sin0 
x  y 


~  w. 


(2.5) 


The  moment  equilibrium  equation  requires  that 


i) 


ds 


+  )  EM  =  0  =  M  +  dM-M-(V  +  dV)  ds  +  w  ds  •  —  sin0 

C  X  t— 


.  ds  Q 
-  W  ds  -r—  COS0 

y  2 


Neglecting  second  order  quantities  gives 


dM 

ds 


V. 


(2.6) 


The  deformed  shape  of  the  rod  has  been  expressed  in  terms  of 
the  function  0(s).  It  is  sometimes  more  convenient  to  express  the 
shape  as  y  =  k(x). 
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These  coordinates  are  related  by  the  following  equations 


x  =  x  +  /  cos0  ds 

o  ; 
s 

o 


y  =  y  +  /  sinfl  ds 

o  J 
s 

o 


dy 

0  =  arc tan  ( ^  ) 
dx 


/ 


s  +  /  /  1 

o  J 

X 

o 


dx 


where  (x  ,  y  )  are  the  coordinates  of  the  rod  at  s  =  s 
o’  J o'  o 


(2.7) 


2. 3  The  Problem  of  Unknown  Deformed  Shape 

In  this  problem  the  free  shape  of  the  rod  <J)(s)  and  the 

distributed  loading  w  and  w  or  w  and  w  and  any  concentrated  loads 

°  x  y  n  t 

are  specified.  The  deformed  shape  of  the  rod  0(s)  must  then  be 
determined.  Integrating  equation  (2.2)  and  equation  (2.3)  gives 


s 

Tsin0  -  Vcos0  -  /  w  ds  =  c.  (2.8) 

J  y  1 

o 

s 

Tcos0  +  Vsin0  +  /  w  ds  =  c„  (2.9) 

o 


If  the  boundary  conditions  T  =  T  ,  V  =  Vq  and  0  =  0  at  s  =  0  are 


» 
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substituted  into  equations  (2.8)  and  (2.9)  then 


c-  =  T  ■  sing  -  V  cos$ 
loo 

(2.10) 

Cr,  =  T  cos3  +  V  sin3 
2  o  o 


Multiplying  equation  (2.8)  by  -  cos6  and  equation  (2.9)  by  sin0  and 
adding  gives 


s 

V  +  sin9  / 
o 


w  ds  +  cosQ 
x 


/ 


o 


w  ds  =  Cr,  sin0 

y  2 


cos0 


Substituting  for  V  from  equation  (2.6)  and  for  M  from  equation  (2.1) 
then  gives 


d_ 

ds 


EI  (~-  -  ) 

ds  ds 


s  s 

+  sin0  I  w  ds  +  cos0  f  w  ds 
J  x  y 


=  sin0  -  cos0 


(2.11) 


Equation  (2.11)  is  a  second  order  nonlinear  ordinary  differential 
equation  for  the  deformed  shape  0(s)  subject  to  boundary  conditions 
for  a  particular  problem. 

If  the  normal  and  tangential  loads  and  w^_  are  specified 
then  equation  (2.11)  becomes  an  integral  differential  equation  for 


0(s)  given  by 


t  »  * 

* 

*  X  • 


d_ 

ds 


EI(di  _  di) 

Ms  ds' 


+  sin0  /  (w  cos0  +  v;  sin0)  ds 
t  n 


m 

+  cos0  I  (w  cos0  -  w  sin0) 
J  n  t 


ds 


=  sin0  -  cos0 


(2.12) 


A  general  numerical  method  is  presented  in  Chapter  III  for 
solving  this  type  of  problem. 

2 . 4  The  Problem  of  Unknown  Free  Shape 

In  this  problem  the  deformed  shape  0(s)  and  the  loading  are 
specified.  The  free  shape  cf>(s)  is  to  be  determined.  Rearranging 
equation  (2.11)  gives 


d 

EI  4L 

_  d_ 

r»  m 

EI^ 

ds 

ds 

»  a 

ds 

ds 

s 

sin0  / 
o 


w  ds 
x 


s 

“  COS0  / 
o 


w  ds  +  c„  sin0 

y  2 


COS0 


(2.13) 


This  is  a  linear  second  order  ordinary  differential  equation 

for  <Ks).  Mitchell  [17]  gives  the  general  solution  of  this  equation 

for  a  cantilever  beam.  If  w  and  are  specified  then  equation  (2.12) 

n  t 

again  gives  a  linear  ordinary  second  order  differential  equation  for 
(p  (s ) .  This  problem  has  been  included  for  completeness  and  will  not  be 


considered  further. 


• 
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2. 5  The  Problem  of  Unknown  Loading 

In  this  problem  the  free  and  deformed  shapes  of  the  rod  are 
specified  and  the  required  loading  must  be  determined.  Substituting 
for  V  in  equation  (2.4)  and  (2.5)  gives 


d0  _  df_ 

ds  "  ds2 


«>£  -  £> 


w  =  0 
n 


(2.14) 


dT  +  _d_ 
ds  ds 


d0 

ds 


+  w 

t 


0 


(2.15) 


These  two  differential  equations  involve  the  unknown  functions 

T(s),  w  fs)  and  w  (s).  Some  additional  relation  between  these  functions 
t  n 

is  required  for  solution  of  this  problem.  For  example  if  the  rod  is 
loaded  by  a  normal  load  wn(s)  only,  then  wt  =  0.  Equation  (2.15)  is 
then  a  first  order  ordinary  differential  equation  which  can  be  solved 
for  the  tension  T(s).  Equation  (2.14)  can  then  be  solved  for  the 
required  normal  load.  If  the  rod  is  loaded  only  by  a  tangential 
distributed  load  then  equations  (2.14)  and  (2.15)  can  be  solved  algebra¬ 
ically  for  the  tension  T(s)  and  the  unknown  tangential  loading  w  . 

Another  more  practical  case  is  that  of  a  rod  loaded  by  a  known 
distributed  load  w  in  the  positive  w^  direction  (for  example  the  weight 
of  the  rod).  Proportional  normal  and  tangential  loads  F(s)  and  pF(s) 
are  then  applied  to  the  rod  to  deform  it  to  the  required  shape. 
Substituting  for  this  loading  in  equations  (2.14)  and  (2.15)  and 


eliminating  F(s)  from  the  resulting  equations  leads  to  a  first  order 


ordinary  differential  for  T(s).  The  normal  load  F(s)  can  then  be 
obtained  from  equation  (2.14). 

Chapter  IV  considers  solutions  of  this  type  of  problem  for 
the  case  of  an  initially  straight  rod.  Analytical  solutions  are  given 
for  special  cases  and  a  numerical  integration  is  used  for  more  general 


shapes . 
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CMP  TER  III 

NUMERICAL  SOLUTION  OF  THE  PROBLEM  OF 
UNKNOWN  DEFORMED  SMPE 

3 . 1  Basic  Numerical  Method 

This  chapter  considers  the  problem  of  determining  the  shape 
of  an  initially  straight  rod  loaded  by  distributed  and  concentrated 
loads.  The  numerical  solution  presented  can  be  extended  to  treat 
initially  curved  rods. 

It  was  found  in  section  2.3  that  a  second  order  nonlinear 
differential  equation  must  be  solved  to  obtain  the  deformed  shape  of 
the  rod  G(s).  The  solution  of  this  nonlinear  equation  can  be  avoided 
by  considering  the  rod  as  a  series  of  short  segments.  Within  each 
segment  only  small  relative  rotations  occur  while  the  segment  as  a 
whole  undergoes  large  rotations.  If  the  problem  is  treated  as  an 
initial  value  problem  then  the  overall  shape  can  be  built  up  by  solving 
a  linear  initial  value  problem  in  each  segment.  The  initial  conditions 
used  in  each  segment  can  be  obtained  by  using  geometric  and  force 
compatibility  conditions  at  the  junction  with  the  previous  segment.  The 
actual  boundary  value  problem  is  solved  by  considering  a  series  of 
overall  initial  value  problems  which  converge  to  the  required  boundary 
conditions. 

The  use  of  these  segments  allows  solution  of  problems  involving 
varying  distributed  loads  and  varying  flexural  rigidity,  as  well  as 


- 


. 
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problems  with  several  concentrated  loads. 

3 . 2  Solution  for  the  Rod  Segment 

Since  each  segment  undergoes  only  small  relative  rotations 
the  angle  0  will  be  small  with  respect  to  some  local  x  and  y  axis  as 
shown  in  Figure  5. 


FIGURE  5:  General  Rod  Segment 


Neglecting  the  term  0  as  small  compared  to  unity  for  an 
initially  straight  rod  equation  (2.11)  becomes 


EI  +  sw  6  +  sw  =  T  (0  -  6)  +  V 
2  x  y  o  o 

ds 


(3.1) 


where  w  and  w  are  uniform  distributed  loads  acting  parallel  to  the 
x  y 


' 


. 
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local  x  and  y  axis  and  the  flexural  rigidity  El  is  constant  along  the 

element.  The  constants  c^  and  c^  in  equation  (2.11)  have  been 

evaluated  using  equation  (2.10)  for  the  initial  conditions  T  --  T  , 

V  =  V  and  0  =  (3  at  s  =  0. 
o 

Equation  (3.1)  is  a  linear  ordinary  differential  equation 

which  can  be  solved  for  the  deformed  shape  of  the  segment  0(s)  subject 

to  the  initial  conditions  M  =  M  and  0  =  3  at  s  =  0. 

o 

This  equation  can  be  nondimens ionalized  by  multiplying  by 

2 

6  where  6  is  the  length  of  the  segment  and  then  dividing  by  the 
flexural  rigidity. 

This  gives 


d20 

dp2 


+ 


w  6 
x _ 

El  " 


P  “ 


T  6 
o 

El 


2 


0  = 


2  3 

V  6  w  6 

o_y 


EI 


EI 


T  L 
o 

EI 


6  (3.2) 


where  p  =  qr  . 

o 

Defining  the  quantities 


T 

O 


T 

o 


EI 


v  = 
o 


V 

o 


EI 


X 


w  6 
x 


EI 


w  6 
EI 


(3.3) 


then  equation  (3.2)  can  be  written  as 
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d2e 

dp2' 


+  [XP  -  T  ]  0  =  v  -  T  3  -  Ipp 


o  o 


(3.4) 


For  an  initially  straight  rod  equations  (2.1)  and  (2.6)  become 


M  =  El 


dO 

ds 


V  =  El 


d2e 

ds2 


The  nondimensional  bending  moment  and  shear  force  are  then 


given  by 


y  = 


M6  =  d0 
El  dp’ 


v  = 


V6‘ 


El 


d2e 

dp2 


(3.5) 


The  solution  to  the  homogeneous  part  of  equation  (3.4)  can  be  found  in 
terms  of  Bessel  functions  however  it  is  more  convenient  to  use  a  power 
series  solution. 

Assuming 


0  =  Z  a  p 
k=0  k 


(3.6) 


then  the  initial  conditions  0-3  and 


M  6 
d0  _  o 

dp  El  Uo 


at  p  =  0 


give 
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a  =  3 
o 


a  =  y 

1  o 


(3.7) 


Substituting  equation  (3.6)  into  the  differential  equation  (3.4)  and 
solving  for  the  coefficient  gives 


a0  =  v  /2 
2  o 


a3  =  (To  y0  “ 


ak  = 


x  a,  0  -  x  ai  o 
o  k— 2  A  k-3 

k(k  -  1) 


k  =  4,  5,  6  . 


(3.8) 


Neglecting  the  term  8  as  small  compared  to  unity  gives  the  local  x  and 
y  coordinates  of  the  deformed  shape  from  equation  (2.7)  as 


x  =  s 


s 

y  =  /  6ds 

o 


where  x  =  0  and  y  =  0  at  s  =  0.  The  nondimens ional  coordinates  are  then 


X  = 
6 


Z 

k=0 


ak  P 


k+1 


k  +  1 


(3.9) 


Substituting  for  c 2  in  equation  (2.9)  from  equation  (2.10)  and  applying 
the  small  angle  approximation  gives 
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T  =  T  +  V3“V0-ws 

O  O  X 


where  w  is  constant, 
x 

Nondimensionalizing  this  equation  gives 


T  +  V  3  -  V0 

o  o 


XP 


(3.10) 


If  uniform  normal  and  tangential  distributed  loads  w  and  w 

n  t 

are  applied  to  the  segment  then  equation  (2.12)  can  be  used.  Applying 
the  small  angle  approximation  to  this  equation  gives 


El 


+ 


S  +  0W 
n 


s 

J*  0ds  + 
o 


w  s 
n 


s 

wt  /  0ds 
o 


=  T  (0  -  3)  +  V 
o  o 

where  again  the  rod  is  initially  straight  and  the  flexural  rigidity  El 
has  been  taken  as  constant  in  the  segment. 

Rearranging  this  equation  and  defining 


s 

0s  =  /  3ds 

o 


El  +  W  (1  +  00)  S  +  w  (0  -  0)  s  =  T  (0  -  3)  +  V 
.  2  n  L  u  u 

ds 


gives 
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The  term  00  can  be  neglected  compared  to  unity.  Tills  gives 

d20  s 

El  — j  +  (w  s  -  T  )0  -  w  /  0ds 
,  z  t  o  t  J 

ds  o 

=  V  -  T  3  -  w  s  (3.11) 

o  o  n  \  / 


Defining 


n 


w  6 
n 


El 


and 


6 


*3 

EI~ 


then  equation  (3.11)  can  be  nondimens ionalized  in  the  same  manner  as 
equation  (3.2)  to  give 

d20  p 

— j  +  (£p  -  x  )  0  -  £  /  0dp 

dP  o 

=  v  -  i  g  -  np  (3.12) 

o  o 


Using  the  power  series  solution  given  by  equation  (3.6)  then  and 
are  given  by  equation  (3.7). 

Substituting  this  power  series  into  (3.12)  yields 


a0  =  v  /2 
2  o 


a.  =  (t  u  -  n)/6 
3  00 


ak 


t  a.  «  -  £  a,  -  (1  -  1/k) 
o  k-2 _ k-3 _ 

k(k  -  1) 


(3.13) 


k  4 )  5 )  6  ... 


■.  - 


. 


* 
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Substituting  for  wx  and  c 2  in  equation  (2.9)  and  using  the 
small  angle  approximation  gives 

s 

T  =  T  +  V  3  -  V0  -  v*  s  ~  w  /  0ds 
o  o  t  n  J 

o 

The  nondimensional  tension  is  then  given  by 

t=t  +  v  6  -  v0  -  £p  -  ny/6  •  (3.14) 

o  o  J 

3 . 3  Equations  for  the  Numerical  Method 

The  rod  is  divided  into  N  segments  each  of  length  6,  with 
the  segments  numbered  from  1  at  the  starting  point  to  N.  The  scheme 
for  assembling  the  segments  is  shown  in  Figure  6 


FIGURE  6:  Scheme  for  Assembling  Segments 


. 


. 
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til 

where  x.  and  y.  are  the  local  coordinate  axis  for  the  j  segment  and 
1  1 

X  and  Y  are  the  global  set  of  axis.  3.  and  cf> .  are  respectively  the 

J  1 

angles  from  the  local  x^  axis  to  the  tangents  at  the  beginning  and 
til 

end  of  the  j  segment.  The  angle  y^.  is  measured  from  the  global  X 
axis  to  the  tangent  at  the  end  of  the  j ^  segment. 


The  bending  moments  at  the  beginning  and  end  of  the  j 


.th 


segment  are  respectively  and  ,  the  tensions  C_.  and  and  the 

shear  forces  F.  and  G.  as  shown  in  Figure  6. 

J  J 

Global  nondimensional  quantities  can  be  defined  as 


A* 

J 


&  - 


A.L 

El* 


B.L 

B*  = 
j  El* 


C.L 

C*  = 


2 


El* 


D.L^ 
D*  =  — 


El* 


(3.15) 


F  L 


G.V 


F*  =  ~J- 


EI* 


G*  =  -1 


El* 


where  L  is  the  length  of  the  rod  and  El*  is  a  representative  flexural 
rigidity  for  the  rod. 

Nondimensional  quantities  based  on  the  segment  length  and 
_  ,  .  th 

flexural  rigidity  of  the  j  segment  are 


a . 
J 


A. 6 

El. 

J 


c .  = 


C.6‘ 

El. 

J 


f  .  = 


F.6‘ 

1 


j  El. 


b.  = 


B.6 

_JL_ 

El. 

J 


d.  = 


D.6 
_J _ 

El . 
J 


8j 


G  .6 
_JL_ 

El. 

J 


(3.16) 


. 
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In  the  numerical  method  it  is  necessary  to  convert  from 
global  nondimensional  quantities  to  segment  nondimensional  quantities 
at  the  beginning  of  the  j ^  segment.  Using  equations  (3.15)  and 
(3.16)  then 


a.  =  A*/(r.  N) 
3  3  3 


c . 
J 


=  C*/(r.  N  ) 
3  3 


(3.17) 


f.  =  F*/(r.  N2) 
J  J  J 


where 


r  . 
3 


El. 
_ 1 

El* 


(3.18) 


and  N  is  the  number  of  segments. 

th 

The  inverse  conversion  is  used  at  the  end  of  the  j  ‘  segment 

given  by 


B*  =  b. 
3  3 


x.  N 
J 


D*  =  d.  r  N2 
j  J  3 


2 

p*  =  o  r  N 
j  j 


(3.19) 
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Before  the  series  solution  for  0 (p )  can  be  evaluated  in  the 

j  segment,  the  angle  3.  which  determines  the  position  of  the  x.  and 

3  3 

y.  axis  must  be  chosen.  The  simplest  choice  is  to  let  3.  be  zero, 

J  J 

then  the  local  axis  will  be  tangent  to  the  beginning  of  the  segment. 

If  the  x^.  axis  was  placed  through  the  end  points  of  the  segment,  then 

the  angle  0(s)  measured  with  respect  to  this  axis  intuitively  would 

be  smaller  at  points  along  the  segment  than  the  angle  taken  with 

respect  to  the  tangential  axis.  In  deriving  the  solution  for  the 

2 

general  segment  in  section  3.2  the  term  0  was  neglected  as  being 

small  compared  to  unity.  It  is  desirable  to  make  this  error  term  as 

2 

small  as  possible  in  the  segment.  The  average  value  of  0  in  the 
segment  is  given  by 


2 

(0  ) 


AVE 


=  /  e2  dp 


(3.20) 


o 


From  equations  (3.6),  (3.7)  and  (3.8)  for  the  general  rod  segment 


v  p  00  k 

0  =  3+  up  +  ~~ —  +  Z  a  p 
°  2  k=3  k 


Approximating  0  by  the  first  three  terms  of  this  series 


gives 


V  ^ 

6  =  3  +  p  p  +  p 

o  l 


(3.21) 


* 


' 


- 


■ 
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and  from  equation  (3.20)  for  this  linear  curvature  approximation  of  the 
shape 


? 

(0  ) 


AVE 


y  v 

+  28  (~°  +  + 


+ 


Minimizing  this  with  respect  to  3  gives 
3y  +  v 

6  =  -  (-- °---°-  )  •  (3.22) 


From  equation  (3.9)  the  local  y-coordinate  of  the  end  point  of  the 
segment  with  the  shape  given  by  equation  (3.21)  is 


y  v 

y  -  6(8  +f  +  f) 


which  equals  zero  when  3  is  given  by  equation  (3.22).  Thus  with  this 

choice  of  3  the  local  x  axis  passes  through  the  end  points  of  the 

linear  curvature  approximation  of  the  deformed  shape  of  the  segment. 

If  only  the  first  two  terms  of  the  series  are  used  then  this 

gives  a  circular  arc  approximation  for  the  shape  of  the  segment. 

2 

Minimizing  0  gives  3  =  -  yQ/2  which  is  equivalent  to  positioning 

the  local  axis  through  the  end  points  of  the  circular  arc  approximation 
of  the  deformed  segment. 

If  only  the  first  term  of  the  series  is  used  then  0=3  and 
o 

0  is  minimized  when  3=0.  The  x-axis  is  then  tangent  to  the 

AVE 


* 
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beginning  of  the  segment.  Using  equation  (3.22)  the  angle 

determines  the  orientation  of  the  local  x.  and  y.  axis  can 

1  3 

as 


3 .  =  -  (3a .  +  f . )/6 
J  J  J 


Evaluating  equations  (3.5)  and  (3.6)  at  the  end  of  the  j 
p  =  1  gives 


b. 

3 


00 


z 

k=l 


E  k(k-l)  a 
k=2 


The  local  coordinates  of  the  end  of  the  segment  are 


and 

\ 


h . 

j 


z 

k=0 


k+l 


3  .  which 
3 

be  chosen 


(3.23) 

segment  where 


(3.24) 


(3.25) 


(3.26) 
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If  the  segment  is  loaded  by  constant  nondimensional  distributed 
loads  Xj  and  ipj  parallel  to  the  x^  and  y.  axis  then  from  equations  (3.7) 
and  (3.8) 


a  =  3  a,  =  a .  a„=f./2 

o  1  J  2  j 


a 


3 


<10/6 

3 


(3.27) 


<cj  ak-2  "  Xj  \-3)/[k(k'1)] 


k  4 j 5 j 6  ... 


and  from  equation  (3.10) 


d. 

J 


c  .  +  f .3  . 
J  3  3 


(3.28) 


If  constant  nondimensional  distributed  loads  and  £  act 
normal  and  tangential  to  the  j  segment  then  from  equations  (3.7)  and 
(3.13) 


a 

o 


=  3 


a  =  (c .  a  -  n  )/6 
J  J  J 

a,  =  [c.  a.  0  -  K.  a  (1-1/k) ]/ [k(k-l) ] 
k  j  k-2  J  k-j 


(3.29) 


k  =  4,5,6 


•  •  • 
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and  from  equation  (3.14) 


d.  =  c.  -  f.3.  -  g.<p.  -  K.  -  n.h. 
3  3  3  3  3  3  3  3  3 


(3.30) 


In  some  problems  both  types  of  constant  distributed  loads  may  be 
applied  to  the  segment.  Since  both  element  solutions  satisfy  linear 
ordinary  differential  equations  the  superposition  principle  can  be 


used  to  give  the  values  of  for  the  combined  loading  of  the  j 


.  th 


segment. 


ttl 

The  angle  y  at  the  end  of  the  j  segment  from  Figure  6 


is  given  by 


Yj  =  Yj_l  +  b  -  b 


(3.31) 


th 

If  X  ,  Y.  are  the  global  coordinates  of  the  end  of  the  j  segment  then 
3  3 

from  Figure  6. 


X.  =  X._x  +  (K.)|p=1  cosn.^  -  e.)  -  (y.)|p=1  sin(Yj_1  -  gp 

b  =  b-i +  (yj)lp-i  cos(b-i  ‘  V  +  xj'p=i  sin(b-i "  V 


th 

The  nondimensional  global  coordinates  of  the  end  of  the  j  segment  are 


then  given  by 


32. 


x*  =  -1  =  X*_x  +  i  [cos (Yj_1  -  e.)  -  h.sinCY.^  ~  B.)] 

Y. 

Y*  =  -L  =  Y*  +  -  [sin (y  -  B.)  +  h.cos(Y.  .  -  8.)]  (3.32) 

J  *-»  J"'1  J  “-L  J  J  J  “J-  j 


tVi 

Given  the  quantities  B  ^  at  the  end  of  the  j  segment 

then  if  no  concentrated  forces  act  on  the  rod,  the  forces  and  bending 

til 

moment  acting  at  the  start  of  the  (j+1)  segment  are  given  by 


B. 

J 


D. 

3 


or  in  nondimensional  form 


B* 

J 


C*  =  D*  (3.33) 

j+1  j 


G* 

J 
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If  concentrated  forces  or  couples  are  applied  to  the  rod 
then  these  are  assumed  to  act  at  the  junction  of  the  segments.  Equation 
(3.33)  must  then  be  modified  to  account  for  these  loads. 

3 . 4  Numerical  Procedure 

The  usual  nonlinear  bending  problem  solved  is  the  two  point 
boundary  value  problem  in  which  boundary  conditions  are  specified  at 
both  ends  of  the  rod.  An  easier  problem  to  solve  is  the  initial  value 
problem  where  all  boundary  conditions  are  specified  at  one  end  of  the 
rod.  The  numerical  procedure  for  solving  this  problem  is  given  in 
subsection  3.4.1.  The  boundary  value  problem  can  be  considered  as  an 
initial  value  problem  in  which  some  of  the  initial  conditions  are  not 
known.  If  these  initial  conditions  can  be  found  so  that  the  solution 
to  the  initial  value  problem  also  satisfies  the  boundary  conditions 
at  the  end  of  the  rod  then  this  constitutes  the  solution  to  the 
boundary  value  problem.  In  subsection  3.4.2  an  iterative  method  is 
presented  which  can  be  used  to  find  these  unknown  initial  conditions 
and  thus  obtain  the  solution  to  the  boundary  value  problem. 

3.4.1.  Numerical  Solution  of  the  Initial  Value  Problem 

If  the  bending  moment,  shear  force,  tension,  global  coordinates 
and  slope  of  the  beginning  of  the  rod  are  known  then  the  solution  for 
the  total  rod  can  be  built  up  by  solving  each  segment  successively 


34. 


starting  with  segment  1.  Figure  7  shows  a  general  rod  with  the  six 
boundary  quantities  at  each  end. 


FIGURE  7 :  Boundary  Quantities  of  Rod 


The  nondimens ional  coordinates  at  the  beginning  of  the  rod 
are  given  by 

X  Y 

X*  =  ^  and  Y*  =  ~ •  (3.34) 

o  L  o  L 

where  L  is  the  length  of  the  rod  and  the  slope  is  given  the  angle  y  . 

In  terms  of  the  notation  for  the  first  segment  the  global  nondimensional 
moment,  shear  force  and  tension  at  the  beginning  of  the  rod  are  A*,  F*, 
C*.  Given  these  six  quantities  then  the  solution  for  the  first  segment 


can  be  obtained. 
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For  the  j  element  given  the  initial  quantities  A*.  F*  C* 

J  J  J 

Y  *  _->  s  X*  -  ,  Y*  and  the  flexural  rigidity  ratio  r.  then  the  following 
procedure  can  be  used  to  obtain  the  quantities  at  the  beginning  of  the 
next  segment. 

1,  The  segment  nondimensional  moment  a.,  shear  force  f.,  and 

J  J 

tension  c^  are  calculated  from  equation  (3.17). 

2.  The  angle  3.  which  determines  the  local  x.  and  y.  axis  is 

J  J  J 

chosen  using  equation  (3.23). 


3.  The  series  solution  given  in  equations  (3.24)  and  (3.26)  is 

used  to  calculate  the  angle  <j>  taken  with  respect  to  the 

local  x_.  axis,  the  segment  nondimensional  moment  b  ,  shear 

ttl 

force  g.,  and  local  coordinate  h.  at  the  end  of  the  j 
j  j 

segment.  The  coefficient  a  in  the  series  solution  are  given 

K 

by  equation  (3.27)  if  nondimensional  uniform  distributed  loads 

X.,  .  are  applied  parallel  to  the  x.  and  y.  axis.  If 

J  3  J  i 

uniform  distributed  tangential  and  normal  loads  and  q 
are  applied  to  the  element  then  the  coefficients  are  given 


.  th 


by  equation  (3.29). 

4.  The  segment  nondimensional  tension  at  the  end  of  the  j 

segment  d.  is  then  given  by  equation  (3.28)  if  X-  and 

J  J  J 

are  applied  and  by  equation  (3.30)  if  loads  £  and  n_.  are 
given. 

The  global  nondimensional  coordinates  of  the  end  of  the  j 

segment  X*,  Y*  are  then  calculated  from  equation  (3.32). 

J  J 


th 


5. 


\ 


' 
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6.  The  angle  y^  defining  the  slope  of  the  end  of  the  j ^  segment 
with  respect  to  the  global  X^.  axis  is  calculated  from 
equation  (3.31). 

7.  The  global  nondimensional  moment  B*.  shear  force  G*  and 

J  J 

tension  D*  at  the  end  of  the  j*’*1  segment  are  then  obtained 
from  equations  (3.19). 

8.  Equation  (3.33)  can  then  be  used  to  calculate  the  global 

nondimensional  moment  A^+^,  shear  force  and  tension 

at  the  beginning  of  the  (j+l)^  segment. 

The  quantities  A* . n ,  F* , _ ,  C* , . ,  y.,  X*,  Y*  are  now  known. 

1+1  1+1  1+1  1  1  1 

The  procedure  can  be  repeated  successively  to  obtain  the  end  point 

quantities  for  all  segments.  The  nondimensional  moment  B*,  shear  force 

G*,  tension  D*  coordinates  X*  and  Y*,  and  angle  y.T  at  the  end  of  the 

ir  N  N  N  N 

til 

N  segment  then  become  the  end  conditions  for  the  rod  as  shown  in 
Figure  7. 


3.4.2.  Solution  of  the  Boundary  Value  Problem 


For  the  initial  value  problem  six  quantities  had  to  be 
specified  at  the  start  of  the  rod  to  give  a  complete  solution.  In 
the  two  point  boundary  value  problem  if  (6_M)  quantities  are  specified 
at  the  beginning  of  the  rod  then  M  quantities  (q^,  ^2*  *  °  *  9^)  must 
be  specified  at  the  end  of  the  rod.  Let  the  specified  values  of  these 
quantities  be  q*,  q*,  ...  q*.  This  means  there  are  M  unknown  quantities, 


o  •  © 
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Z-j^j  ^*2 *  *  *  *  beginning  of  the  rod  which  must  be  determined. 

In  vector  form  these  quantities  are  given  by 


£  = 


(3.35) 


The  vector  c[  can  be  treated  as  a  function  of  Z  since  for  given  values 
of  the  values  of  c[  can  be  obtained  by  solving  the  initial  value 
problem  as  in  section  3.4.1.  The  boundary  value  problem  can  then  be 
solved  by  finding  values  of  Z_  which  satisfy  the  equation 


£(Z)  =  £(Z)  -  =  0  (3.36) 

By  beginning  with  the  end  of  the  rod  which  gives  the  least 
number  of  unknown  initial  conditions  then  at  most  equation  (3.36)  will 
by  a  system  of  three  equations  involving  three  unknowns.  By  tactful 
positioning  of  the  global  axis,  use  of  overall  equilibrium  equations 
and  of  symmetry  of  loading  and  boundary  conditions,  the  number  of 
unknown  initial  quantities  can  be  reduced  to  one  or  two  in  most 
problems . 

This  system  of  equations  can  be  solved  numerically  using 
a  variation  of  Newtons  Method.  Newtons  iteration  is  given  by 


' 
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T-i  j 
1  R 


(3.37) 


where  _p  -  ]3(Z  )  and  where  the  M  x  M  matrix  has  the  general  term 


In  this  problem  the  above  partial  differentiation  cannot  be 
evaluated.  A  backward  difference  formulation  can  be  used  instead.  The 
general  term  of  can  then  be  given  as 


Pk(zb 


Vb1’ 


Z£- 1’  Z£ 


1 1.  z1 


£+1’ 


*  V 


k£ 


y1  -  Z i"1 


(3.38) 


Starting  with  guesses  7?  and  _Z  then  equations  (3.37)  and 
(3.38)  can  be  iterated  until  convergence  is  obtained.  The  one 
dimensional  form  of  this  method  is  usually  called  the  secant  method. 
In  this  work  the  M-dimensional  form  of  this  method  will  be  referred 
to  as  the  M-dimensional  secant  method. 


3. 5  Example  Boundary  Value  Problems 

This  section  presents  four  boundary  value  problems  solved 
using  the  numerical  method  outlined  in  the  previous  section.  The 
first  problem  considered  is  that  of  a  uniform  horizontal  cantilever 
with  a  horizontal  concentrated  load  applied  at  the  free  end.  The 
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numerical  solution  with  various  numbers  of  segments  is  compared  to  the 
known  analytical  solution.  The  problem  of  a  uniform  cantilever  subject 
to  a  uniform  normal  load  is  also  considered.  Again  comparison  is  made 
to  the  analytical  solution.  The  third  problem  is  that  of  a  rod  of  uniform 
weight  per  unit  length  pinned  at  both  ends.  No  analytical  solution  of 
this  problem  has  been  found.  Experimental  verification  of  the  numerical 
solution  of  this  problem  is  given  in  Chapter  V.  The  last  problem 
considered  is  that  of  a  cantilever  of  nonuniform  flexural  stiffness 
subject  to  a  nonuniform  distributed  weight  loading.  The  numerical 
solution  is  compared  to  a  previous  analog  computer  solution.  The  same 
problem  with  an  additional  concentrated  end  load  is  also  considered. 

These  problems  are  solved  using  the  general  computer  program 
given  in  Appendix  I.  Several  lines  of  this  program  must  be  specialized 
to  solve  a  particular  problem.  Appendix  I  gives  the  specialized  groups 
of  lines  used  to  solve  the  above  problems . 

3.5.1.  The  Elastica  Problem 

This  section  considers  the  elastica  problem  of  an  initially 
straight  horizontal  cantilever  loaded  by  a  concentrated  horizontal 
load  P  at  the  free  end  as  shown  in  Figure  8a.  A  two  point  nonlinear 
boundary  value  problem  must  be  solved  to  determined  the  deformed  shape 
of  the  rod. 

This  problem  can  be  solved  analytically  in  terms  of  elliptic 
integrals  (see  Frisch-Fay  [9]).  With  axis  located  as  in  Figure  8a,  the 
coordinates  of  the  deformed  shape  are  given  by 


' 
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FIGURE  8:  The  Elastica 
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X*  =  X/L  =  [2E(<J),p)  -  F(cj>,p)]/K(p)  (3.39) 

Y*  =  Y/L  =  2p  cos<J>/K(p) 

am[pK(p),  p] 

(3.40) 

where  El*  is  the  flexural  rigidity  of  the  rod.  E(cj),p)  and  F(<f>,p)  are 
the  incomplete  elliptic  integrals  of  the  first  and  second  kinds  with 
modulus  p  and  amplitude  cj>.  The  complete  elliptic  integrals  of  the 

TT  TT 

first  and  second  kinds  are  K(p)  =  EC^-,  p)  and  E(p)  =  E(— ,  p). 

The  lengths  d  and  h  as  shown  in  Figure  8a  are  then  given 
from  equation  (2.39)  by 


where  L  is  the  length  of  the  cantilever,  p  =  sin(y/2),  cf> 
and  p  =  s/L.  The  nondimens ional  load  is  given  by 


p*  =  P— ■  =  K2(p) 

El* 


=  2p/K(p) 

p=0 


=  [2E(p)  -  K(p)]/K(p) 

p=l 


(3.41) 


The  numerical  method  given  in  section  3.4  for  the  solution 
of  boundary  value  problems  can  be  used  in  several  ways  to  solve  this 
problem.  With  the  axis  positioned  as  in  Figure  8a  the  X  coordinate 
and  slope  are  the  only  quantities  specified  at  the  clamped  end.  At 


\ 
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the  free  end  of  the  cantilever  the  X  coordinate  and  slope  are  not 
known.  By  starting  at  the  free  end  the  problem  can  then  be  solved 
using  the  two  dimensional  form  of  the  secant  method  given  in  section 
3.4.  If  the  global  coordinate  axis  are  positioned  as  in  Figure  8b  then 
the  slope  at  the  free  end  is  the  only  unknown  quantity.  The  one 
dimensional  secant  method  can  then  be  used  to  solve  the  equation 

p(z)  =  q (Z)  -  q*  =  0  (3.42) 


where  q  is  the  slope  at  the  clamped  end  with  specified  value  q*  -  0  and 

Z  =  y  the  unknown  angle  at  the  free  end  of  the  cantilever.  In  terms 

of  the  notation  given  in  Figure  7  for  the  end  conditions  of  the  general 

rod  then  the  known  initial  conditions  are  X*  =  0,  Y*  =  0,  the  moment 

o  o 

A*  =  0,  the  shear  force  F*  =  -P*siny  and  tension  =  -  P*cosy.  The 

unknown  initial  angle  is  y  =  y  and  the  angle  at  the  clamped  end  of 

o 

the  rod  is  yXT,  with  specified  value  y  =  0. 

N  N 

The  numerical  results  for  this  problem  were  obtained  using 

the  general  computer  program  in  Appendix  I.  For  the  case  of  a 

cantilever  of  uniform  cross  section  the  flexural  rigidity  ratio  for 
i  . th  ^  . 

the  j  segment  is 


r . 
J 


EX. 

_ 1  = 

El* 


1. 


The  analytical  solution  vas  evaluated  using  the  elliptic  integral  tables 
given  in  reference  [23]. 
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In  Table  1  the  values  of  y,  d*  and  h*  given  by  these  solutions 

are  compared.  The  subscripts  N  and  A  refer  respectively  to  the  numerical 

solution  using  N  segments  and  to  the  analytical  elliptic  integral  solution. 

% 

In  the  elliptic  integral  solution  the  angle  y  is  first  specified  and 
then  P*  is  calculated.  This  value  of  P*  is  then  used  in  the  numerical 
method. 

The  program  was  run  with  the  local  segment  axis  positioned 
through  the  end  points  of  the  linear  curvature  approximation  to  the 
segment  of  the  rod,  through  the  circular  arc  approximation,  and 
tangential  to  the  beginning  of  the  segment.  Using  the  linear  curvature 
approximation  to  position  the  local  axis  tended  to  give  the  most 
accurate  results.  However  as  the  number  of  segments  is  increased  the 
differences  in  the  solutions  using  these  three  methods  become  insignifi¬ 
cant.  Also  when  the  number  of  segments  is  doubled  the  error  in  the 
numerical  solution  is  approximately  halved.  This  relationship  becomes 
more  exact  as  the  number  of  segments  is  increased.  If  this  relationship 
is  satisfied  exactly  then 


(3.43) 


where  q* represents  one  of  the  quantities  y,  d*  or  h*  and  the  subscripts 
indicate  the  number  of  segments  used  in  the  numerical  solution.  With 
P*  =  4.65056  then  h*Q  =  0.11908  and  h*6Q  =  0.12111,  using  equation  (3.43) 
gives  a  better  estimate  of  h*  as 
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TABLE  1:  Comparison  of  the  Numerical  Method 
using  Different  Numbers  of  Segments 
to  the  Analytical  Elastica  Solution 


p*  = 

2.55406 

3.  =  -  (3a.  +  f. 
3  3  3 

Number 

of 

yn 

d* 

N 

Segments 

1 

0.65770 

0.41022 

0.91233 

5 

0.55716 

0.34392 

0.92412 

10 

0.54064 

0.33403 

0.92812 

20 

0.53211 

0.32896 

0.93026 

40 

0.52784 

0.32642 

0.93134 

80 

0.52572 

0.32516 

0.93189 

160 

0.52466 

0.32452 

0.93216 

2q160 q80 

0.52360 

0.32388 

0.93243 

Analytical 

Result 

Y .  =  0.52360 

A 

d*  =  0.32390 

A 

h*  =  0.93243 

A 

p*  = 

4.65056 

1 

2.27013 

1.00447 

-0.01135 

5 

2.15960 

0.81142 

0.07114 

10 

2.13035 

0.80572 

0.09312 

20 

2.11305 

0.80409 

0.10733 

40 

2.10387 

0.80354 

0.11507 

80 

2.09916 

0.80334 

0.11908 

160 

2.09679 

0.80325 

0.12111 

2ql60~q80 

2.09442 

0.80316 

0.12314 

Analytical 

Result 

y  =  2.09440 

A 

d*  =  0.80317 

A 

h*  =  0.12316 

A 

)/6 
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h*  =  2(. 12111)  -  .11908  =  .12314 

This  result  is  accurate  to  four  significant  figures  while  h*  and 

80 

hJ^Q  are  only  correct  to  two.  The  values  of  y ,  d*  and  h*  given  by 
equation  (3.43)  with  N  =  80  are  shown  in  Table  1.  This  type  of 
convergence  was  obtained  in  all  problems  involving  uniform  rods  with 
concentrated  and  uniformly  distributed  loads. 

Table  2  compares  the  values  of  the  nondimens ional  X  and  Y 
coordinates  given  by  the  numerical  method  (N  =  160)  to  the  analytic 
elastica  solution  at  five  points  along  the  cantilever.  For  P*  =  2.55406 
the  rod  is  approximately  parallel  to  the  X  axis.  In  this  case  the 
relative  error  in  the  values  of  X*  decrease  along  the  rod  while  the 
errors  in  Y*  increase.  For  the  case  of  P*  =  4.05060  the  rod  is 
approximately  parallel  to  the  Y  axis  and  the  errors  in  Y*  decrease 
along  the  rod  while  the  errors  in  X*  increase. 
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TABLE  2:  Comparison  of  the  Numerical  Method 

to  the  Analytic  Elastica 

Solution 

at  Various  Points  along 

the  Rod 

P*  =  2.55406 

N  =  160  3.  =  -  (3a. 

J  3 

+  f  )/6 

s 

Numerical 

Method 

Analytic  Solution 

P  =  L 

X* 

160 

Y* 

160 

Y* 

A 

0.2 

0.17394 

0.09871 

0.17404 

0.09853 

0.4 

0.35256 

0.18857 

0.35275 

0.18822 

0.6 

0.53888 

0.26095 

0.53913 

0.26046 

0.  8 

0.73306 

0.30812 

0.73325 

0.30753 

1.0 

0.93216 

0.32453 

0.93243 

0.32390 

P*  =  4.65056 

0.2 

-0.09569 

0.17555 

-0.09526 

0.17579 

0.4 

-0.16122 

0.36396 

-0.16030 

0.36437 

0.6 

-0.16026 

-0.56281 

0.15877 

0.56281 

0.8 

-0.06132 

0.73267 

-0.05937 

0,73281 

1.0 

0. 12111 

0.80325 

0.12316 

0.80317 

‘ 
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3.5.2.  The  Cantilever  with  Uniform  Normal  Load 

This  problem  is  presented  to  compare  the  solution  using  the 
numerical  method  to  one  of  the  few  analytic  solutions  obtained  for 
rods  loaded  by  distributed  loads.  The  problem  of  an  initially 
straight  cantilever  loaded  by  a  uniform  normal  load  per  unit  length  P 
has  been  solved  by  Mitchell  [17]. 


FIGURE  9:  Cantilever  with  Normal  Load 


In  terms  of  the  polar  coordinates  shown  in  Figure  9  the  deformed  shape 
of  the  cantilever  is  given  by 


0 


4  .  -1 

‘P*(r*)3 ' 

1  .  -1 

**  m 

P*(r*)3 

i sin 

L  8 

8 

(3.44) 


where 
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p* 


PL 

El* 


and 


The  nondiraensional  load  P*  and  the  radius  r*  are  related  by  the  equations 


F(fi,k)  =73  (P*) 


1/3 


(3.45) 


where  F(fi,k)  is  the  incomplete  elliptic  integral  of  the  first  kind  and 
sin  ^(k)  =  15°.  Also 


Q 


-1 

=  cos 


1  -  (/3  +  l)(P*)2/3(r*)2/2 
1  +  (/3  -  l)(P*)2/3(rT*)2/2 


(3.46) 


These  equations  can  be  solved  by  choosing  fi,  solving  equation  (3.45) 
for  P*  and  then  using  equation  (3.46)  to  determine  r*. 

L 

Then 


Y  =  0  | 


4 


’p*(r*) 

JLi 


p=0 


3  Sln 


3. 


(3.47) 


h*  =  I  =  rL  COS  ®L 


d*  "  l  "  vt  Sln  ®L 


where 


49. 


sm 


-1  P*(rL} 


8 


As  in  the  elastica  problem  this  problem  can  be  solved  using 

the  one  dimensional  secant  method  starting  from  the  free  end.  The 

known  initial  quantities  are  =  0,  =  0,  A^  =  0,  C  =  0,  F  =  0. 

The  angle  y  =  y  must  then  be  found  such  that  the  angle  y  at  the 

o  N 

clamped  end  of  the  cantilever  is  given  by  y  =0. 

N 

The  nondimens ional  normal  load  on  the  j  segment  of  the 

rod  is 

3 

Pa  3 

n.  ~  -  PVV/) 

where  the  flexural  rigidity  ratio  r_.  =  1  for  a  cantilever  of  uniform 
cross  section. 

The  values  of  y,  d*  and  h*  from  the  numerical  method  using 
the  linear  curvature  approximation  to  position  the  segment  axis  are 
compared  to  the  analytical  results  in  Table  3.  This  method  tended  to 
give  better  results  than  using  a  circular  arc  approximation  to  position 
the  local  axis  or  using  tangential  segment  axis.  As  in  the  elastica 
problem  as  the  number  of  segments  is  increased  the  differences  in  the 
results  using  these  three  methods  become  insignificant.  Again  for 
large  numbers  of  segments  the  errors  tend  to  be  halved  when  the  number 


of  segments  is  doubled. 


V  . 
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TABLE  3:  Comparison  of  the  Numerical  Method 
with  the  Analytical  Solution  for  a 
Cantilever  with  Uniform  Normal  Load 


P*  =  3.29814 


Number 

of 

Segments  i 

yn 

d* 

N 

h§ 

1 

0.54969 

0.405244 

0.92448 

2 

0.54852 

0.397524 

0.91181 

3 

0.54788 

0.39567 

0.90829 

4 

0.54749 

0.39510 

0.90720 

5 

0.54720 

0.39464 

0.90676 

10 

0.54641 

0.39385 

0.90637 

20 

0.54590 

0.39340 

0.90638 

40 

0.54560 

0.39317 

0.90647 

Analytical 

Result 

y .  =  0.54530 

1 A 

d*  =  0.39292 

A 

h*  =  0.90657 

P*  =  14.32502 


1 

2.38750 

1.11968 

-0.32026 

3 

2.24554 

0.82766 

-0.09431 

5 

2.20156 

0.81722 

-0.07390 

10 

2.15238 

0.81810 

-0.04230 

20 

2.12397 

0.82041 

-0.02205 

40 

2.10926 

0.82188 

-0.01119 

80 

2.10184 

0.82267 

-0.00563 

160 

2.09811 

0.82308 

-0.00282 

Analytical 

Result 

Y  =  2.09440 

A 

d*  =  0.82350 

A 

h*  =  0.00000 

A 

4 
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3.5.3.  The  Stiffened  Catenary 

The  deformed  shape  of  a  rod  of  uniform  weight  per  unit 
length  W  when  pinned  at  both  ends  has  not  been  determined  analytically. 
Given  the  coordinates  of  the  ends  of  the  rod  and  that  the  bending 
moment  at  each  end  is  zero  then  the  three  dimensional  secant  method 
could  be  used  to  determine  the  three  remaining  initial  conditions. 

For  the  symmetric  case  shown  in  Figure  10a  the  vertical  component  of 
the  reaction  at  each  end  of  the  rod  is  WL  where  L  is  the  half  length 
of  the  rod.  Using  this  fact  then  the  problem  could  be  solved  using 
the  two  dimension  secant  method. 

In  Chapter  V  this  problem  is  treated  experimentally  with 
the  ends  of  the  rod  supported  by  strings  at  some  angle  a  with  respect 
to  X  axis  as  shown  in  Figure  10b.  If  the  global  axis  are  positioned 
as  in  Figure  11  then  the  solution  can  be  found  using  a  one  dimensional 
secant  method  starting  at  the  pinned  end  of  the  rod. 

The  known  initial  conditions  are  X*  =  0,  Y*  =  0,  and  the 

o  o 

moment  A*  =  0.  From  symmetry  the  shear  force  at  the  mid  point  of  the 
rod  is  zero.  For  equilibrium  the  vertical  component  of  F  is  F^  =  WL 
and  the  horizontal  component  is  F  =  -  WL  ctn  a.  The  tension  and 
shear  force  at  the  beginning  of  the  rod  are  then 

C1  =  WL(siny  +  etna  cosy) 

F1  =  WL(cosa  -  etna  siny) 


* 
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Y  A 


*x 


(a) 


(b) 

FIGURE  10:  The  Stiffened  Catenary 


FIGURE  11:  Position  of  Global  Axis 
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In  global  nondimensional  form 


C-^*  =  W*(siny  +  etna  cosy) 


F^*  =  W*(cosy  -  etna  siny) 


where  W*  is  the  nondimensional  weight  of  half  the  rod 


W*  = 


WL 

El* 


where  the  reference  flexural  rigidity  El*  is  taken  as  the  flexural 
rigidity  of  the  rod. 

The  angle  y^  =  -  y  is  the  only  unknown  initial  quantity. 

A  one  dimensional  secant  method  can  be  used  to  find  y  so  that  one 

o 

of  the  known  boundary  conditions  at  the  end  of  the  rod,  say  y^  =  0 


is  satisfied. 


th 


The  distributed  loads  used  in  the  j  element  are 


w  =  -  W  sin(y  -  6.) 
x  J-l  J 


w  =  W  cos  (y .  ..  -  3  . ) 
y  j-1  3 


The  nondimensional  loads  are  then 


• 

* 

V 

, 
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*3 

w  o 

X.  =  ~~ —  =  "  W*  sin (y  -  3.)/(r.N3) 
J  EIj  3  J 

'  *3 

w  o 

^  =  -^T—  =  W*  cos  (y  -  3.)/(r.N3) 

J  Eij  J  3 


For  a  uniform  rod 


The  numerical  results  were  obtained  using  the  general  program  in 
Appendix  I  specialized  for  this  problem.  Figure  12  shows  the  shape  of 
the  rod  for  various  nondimensional  weights  with  a  =  arctan  2.  Keep¬ 
ing  the  angle  a  constant  and  increasing  the  nondimensional  weight 
W*  is  equivalent  to  decreasing  the  stiffness  of  the  rod  while  keeping 
the  actual  weight  of  the  rod  and  the  end  reactions  constant.  As  W*  is 
increased  the  shape  of  the  rod  approaches  the  true  catenary  solution. 

Figure  13  can  be  used  to  determine  the  dimensionless  vertical 
deflection  d*  =  d/L  and  the  horizontal  span  h*  =  h/L  for  rods  over 
the  entire  range  of  possible  nondimensional  weights,  and  for  values 
of  the  angle  a  ranging  from  a  =  0  to  a  =  ir/2.  The  case  of  a  rod 
supported  only  by  vertical  end  reactions  occurs  when  a  =  tt/2.  As  a 
is  decreased  the  horizontal  component  of  the  end  reaction  increases. 

This  horizontal  component  approaches  infinity  as  a  approaches  zero. 


-  .  *  • 


, 


. 


* 
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FIGURE  12:  Shapes  of  Stiffened  Catenary  with  Varying 

Flexural  Stiffness 
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FIGURE  13:  Vertical  Deflection  and  Horizontal  Span  of 

Stiffened  Catenary 
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3.5.4.  The  Tapered  Heavy  Cantilever  with  Concentrated  Load 

In  this  problem  a  trapezoidal  strip  of  uniform  density  is 

deformed  by  its  weight  and  a  concentrated  load  P  at  its  free  end. 

The  load  P  will  be  treated  as  a  follower  force  acting  at  a  specified 

angle  a  with  respect  to  the  cantilever  as  shown  in  Figure  14a.  It 

could  also  be  considered  as  acting  at  a  specified  angle  to  some 

fixed  axis  (for  example  the  angle  a  +  y  could  be  kept  constant). 

This  problem  can  be  solved  numerically  using  a  one  dimensional 

secant  method  starting  at  the  free  end.  The  known  initial  conditions 

are  X  =  0,  Y  =  0,  A*  =  0,  C*  =  P*  cosa  and  F*  =  -P*  sina,  where 
o  o  1  1  I 

2 

P*  =  PL  / El*.  The  reference  flexural  rigidity  El*  is  taken  as  the 

flexural  rigidity  at  the  base  of  the  cantilever  where  s  =  L. 

The  top  view  of  the  cantilever  is  shown  in  Figure  14b.  The 

ratio  of  the  widths  is  defined  as  t  =  b  /b  .  If  t  =  0  then  the 

o  L 

cantilever  is  triangular.  The  special  case  of  a  uniform  cantilever 
is  given  by  t  =  1. 

The  width  of  the  strip  at  any  point  is  given  by 
b(p)  =  bL[t  +  p(l-t)]  (3.48) 

where  p  =  s/L. 

For  a  strip  of  constant  thickness  the  flexural  rigidity  at 
any  point  is  proportional  to  the  width  of  the  strip.  This  gives  the 


flexural  rigidity  ratio 


• 

, 

■v  ** 

/ 
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TOP  VIEW 


FIGURE  14:  The  Tapered  Cantilever 


* 
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r(p) 


=  EI(P) 

El* 


t  +  p(l-t) 


(3.49) 


The  weight  per  unit  length  of  the  strip  W  will  also  be  proportional  to 
the  width  b(s)  giving 


W(s) 

W_ 


t  +  p  (1-t) 


(3.50) 


where  W  is  the  weight  per  unit  length  of  the  strip  at  p  =  1. 

In  the  numerical  procedure  segments  with  constant  flexural 

rigidity  and  uniform  loads  must  be  used.  The  shape  is  approximated 

by  constant  width  segments  with  the  width  taken  at  the  mid  point  of 

the  segment.  Starting  at  the  free  end  of  the  strip  with  segment  1, 

t  h 

the  mid  point  of  the  j  segment  is  given  by 


Pj  =  (j-.5)/N 


The  flexural  rigidity  ratio  for  the  j 


th 


element  is  then 


r . 
J 


El* 


r  (p  j  ) 


(3.51) 


(3.52) 


The  distributed  loads  in  the  segment  are 


w  =  -  W(p  .)  sin  (y  ,  -  3  . ) 
x  J  J-l  J 


. 
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wy  =  w^pj^  cos(^Yj-i  “ 


In  nondimens ional  form  after  simplification  using  equations 
(3.50),  (3.51)  and  (3.52)  then 


sin <Yj_i  -  Bp/N3 


*3 


w  6 ' 

_X _ 

El. 


W*  cos  (y  .  -  3  .  )/N~ 

f  J"1  3 


where  W*  -  W  L  /El*  is  the  nondimensional  weight  per  unit  length  at 

i-!  i-J 

the  base  of  the  cantilever. 

The  general  program  specialized  for  this  problem  as  given 
in  Appendix  I  was  run  for  the  special  case  of  a  triangular  cantilever 
(t  =  0)  with  a  variable  load  P  applied  perpendicular  to  the  free  end 
(a  =  7t/2).  The  results  with  W*  =  10  are  shown  in  Figure  15.  When  P* 

L 

is  zero  the  results  can  be  compared  to  the  analog  computer  solution 

of  Lippman,  Mahrenholtz  and  Johnson  [16].  Their  results  are  tabulated 

for  specified  values  of  the  angle  y.  Given  this  initial  angle  Yq  =  y 

then  the  loading  parameter  W*  can  be  treated  as  the  unknown  quantity 

.Li 

which  must  be  found  using  the  one  dimensional  secant  method  to  give 
a  zero  slope  at  the  base  of  the  cantilever.  Comparison  of  these 
results  is  given  in  Table  4. 


. 


,■  • 
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FIGURE  15:  Solution  for  the  Tapered  Cantilever  with  Varying 

End  Load 
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TABLE  4:  Comparison  of  Numerical  Method 


with  Lippman 

et  al 

[16] 

Y  = 

10° 

Number  of 

Segments 

)r _ 

h*  =  h/L 

d*  =  d/L 

W* 

Li 

10 

0.99033 

0.12987 

3.1081 

20 

0.99025 

0.13033 

3.1549 

40 

0.99023 

0.13045 

3.1669 

80 

0.99022 

0.13047 

3.1700 

Reference [16  J 

0.993 

0.129 

3.18 

Y  = 

30° 

10 

0.91357 

0.38018 

10.060 

20 

0.91288 

0.38130 

10.221 

40 

0.91271 

0.38157 

10.266 

80 

0.91266 

0.38163 

10.278 

Reference [ 16] 

0.899 

0.372 

10.2 

Y  = 

1 

O 

O 

1  ^ 

10 

0.53817 

0.78444 

39.559 

20 

0.53563 

0.78476 

40.354 

40 

0.53511 

0.78468 

40.614 

80 

0.53506 

0.78457 

40.708 

Reference [ 16 ] 

0.539 

0.784 

40.7 

' 


- 
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CHAPTER  IV 

SOLUTION  OF  THE  PROBLEM  OF  UNKNOWN  LOADING 

4 . 1  Basic  Equations 

The  problem  of  determining  the  distributed  loading  necessary 
to  bend  a  rod  to  a  specified  deformed  shape  was  introduced  in  section 
2.5.  The  special  case  of  an  initially  straight  uniform  rod  is 
considered  in  this  chapter.  The  solutions  obtained  could  be  extended 
to  consider  initially  curved  rods. 

For  many  problems  the  loads  are  applied  normal  and  tangential 
to  the  rod.  As  well  one  load  which  often  need  be  considered  is  that 
due  to  gravity.  For  these  cases  the  loads  w^  and  w^  shown  in  Figure 
16  are  given  by 

w  ~  F  +  W  cos0  (4.1) 

n 

w  =  G  -  W  sin9  (4.2) 

where  W  is  the  weight  per  unit  length  and  F  and  G  are  normal  and 
tangential  forces  per  unit  length  applied  to  the  rod. 


. 

' 


■ 
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FIGURE  16 :  Loading  of  Rod 


The  equations  (2.14)  and  (2.15)  using  (4.1)  and  (a. 2)  become 


-  *T  <Li  =  , 


T  — —  El  =  F  +  W  cos0 
ds  ,3 

ds 


(4.3) 


f  +El4f  " 

ds  ds2  ds 


=  -  G  +  W  sin0 


(4.4) 


where  the  flexural  rigidity  El  is  constant.  Equations  (4.3)  and  (4.4) 
can  be  nondimensionalized  to  give 


T 


d0 

dp 


f  +  a)  cos0 


(4.5) 
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dr  d^G  d0 

dp  ,  2  dp 

dp 


-  g  +  a)  sin 6 


where  the  non dimensional  quantities  used  are 


a)  - 


f  = 


WL 

El 

FL 

El 


3 


3 


(4.6) 


(4.7) 


and  in  which  L  is  the  length  of  the  rod.  The  nondimens ional  moment 
and  shear  force  are 


d9 

dp 


and 


(4.8) 


In  equations  (4.5)  and  (4.6)  the  deformed  shape  0(p)  is  known.  These 
equations  contain  three  unknowns  i(p),  f(p)  and  g (p )  and  therefore 
some  additional  relationship  between  these  quantities  must  be  known 
in  order  to  obtain  a  solution. 


4. 2  Proportional  Loading 

While  there  are  numerous  ways  in  which  g  and  f  may  be  related 
and  thereby  allow  solution  of  equations  (4.5)  and  (4.6),  one  possible 
relationship  is 


8  -  yf 


(4.9) 


\ 


' 


1 
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where  y  is  a  constant.  This  can  correspond  to  the  case  of  frictional 
resistance  such  as  would  occur  in  a  bandsaw  blade  with  impending 
sliding  or  during  sliding. 

Substituting  equation  (4.9)  into  equation  (4.6)  and  then 
eliminating  f  from  equation  (4.5)  and  (4.6)  gives 


dx  d0 

--  +  y  —  t 
dp  dp 


d20  dG  ,  .  .  ,  q  ,  d30 

+  co  sin0  +  yco  cos0  +  y  — 


dp 


2  dp 


(4.10) 


dp,' 


The  solution  to  this  first  order  linear  differential  equation  is 


x  =  e 


~y  9 


r  _  ,  f  y0  .  d20  d0  ,  .  .  .  .  yd30.  1  , 

[C  +  I  e  (-  — r  1 - h  co  sin0  +  yco  cos  +  — —)  dp] 

1  J  2  dp  ,  -3 

dp 


dp 


(4.11) 


The  distributed  load  f  is  then  given  from  equation  (4.5)  as 


d©  d  0 

f  =  t  —  -  — r  -  oo  cos0 
dp  ,3 
dp 


(4.12) 


The  constant  C  can  be  evaluated  using 


t  = 


l 

o 


p=o 


Equation  (4.3)  becomes 


T 


e 


y0 

e 


d20  d0  .  n  n  .  yd3©.,. 

(-  — —  h  co  sin0  +  yco  cos0  +  — ^ — ) d p  j 

,  2  dp  .  -J 

dp  dp 


(4.13) 


' 
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The  evaluation  of  equations  (4.11)  or  (4.13)  and  (4.12)  gives  the 
tension  in  the  rod  and  the  normal  force  distribution.  The  moment 
and  shear  force  can  be  evaluated  directly  from  equation  (4.8). 

The  special  case  of  normal  loading  alone  is  given  with 
y  =  0  in  the  previous  solution.  For  the  case  of  only  tangential 
loading  the  normal  load  f  is  zero.  From  equation  (4.5)  the  tension 
in  the  rod  is  given  by 


,  „  .  d~Y.de 

x  =  (co  cos0  +  — y)  — 

,  3  dp 
dp 


Solving  equation  (4.6)  then  gives  the  tangential  load 


g  =  oo  sin0 


d20  d0_  _  dr 

2  dp  dp 
dp 


(4.14) 


(4.15) 


4 . 3  Special  Solutions 

In  this  section  some  special  solutions  of  equations  (4.12) 
and  (4.13),  and  equations  (4.14)  and  (4.15)  are  considered. 


4.3.1.  Circular  Arc 

If  the  rod  is  bent  into  the  circular  arc  then 


0  =  ap 


(4.16) 
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where  a  is  a  constant  which  corresponds  to  the  nondimens ional  moment  ft. 
The  shear  force  is  zero  everywhere. 

For  proportional  loading  equation  (4.13)  can  be  easily 
integrated  to  give 


-yap  .  OJ  r„  .  .  2X  ,  -yap 

t  =  x  e  + - —  L2  ysm  ap  -F  (1-y  )  (e  -  cos  ap)] 

o  .  2. 


a (1+y  ) 


(4.17) 


with  the  corresponding  expression  for  f  from  equation  (4.12)  given  by 


f  =  cxt  -  co  cos  ap 


(4.18) 


If  only  a  tangential  load  is  applied  to  the  rod  then  from  equation 

(4.14) 


co 

t  —  cos  ap 
a 


(4.19) 


and  from  equation  (4.15)  the  tangented  distributed  load  is  given  by 


g  =  2  co  sin  ap 


4.3.2.  Linear  Variation  of  Curvature 


(4.20) 


For  this  case  the  shape  of  the  rod  can  be  given  by 


. 
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ft  ~  ft 

0  =  -  -7T  (1  -  p)  +  T~  +  6  (4 •  21) 

Z  2  o 

where  the  nondimensional  moment  at  p  =  0  is  ft  and  0  is  the  angle. 

o  o 

It  has  been  assumed  in  equation  (4.21)  that  at  p  =  1  the  curvature 
is  zero.  The  nondimensional  shear  force  in  this  case  is  a  constant 
given  by 


v  =  -  ft 


o 


For  proportional  loading  equation  (3.13)  gives 


Pi  y  0  P  o 

t  =  e  ^  (x  e  °  +  /  e^  [ft  (1-p)  +  cosin0  +  ycoccs0]  d0 } 
o  J  o 

o 


(4.22) 


where  0  is  given  by  equation  (4.21).  This  can  be  reduced  to 


t  = 


ft  -yft  p(2-p)/2  ft 

(t  -  — )  e  +  — 

o  y  y 


yft  (1-p)  / 2 
0)  o 

+  -  e 


/2ft~ 


o 


ft  (1-p)  /2 
o '  -yu 

r  /  c  ■  N  t  e  cos  u  , 

[(smc  +  ycosc)  ]  - ryy -  du 

ft  /  2  uU 

o 


ft  (1-p)  /2 
o  -yu  . 

v  r  e  sin  u  ,  ! 

+  (ycosmc  +  cosc)  J  duJ 

ft  /  2  u±/Z 

o 


(4.23) 


ft 


o 


+  0 


where  c  = 


2 


« 


*  < 


. 
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The  remaining  integrals  can  be  evaluated  using  a  power  series 


to  give 


ft  -yft  p (2— p ) / 2 

t  =  (t  - - )  e 

o  y 


yft  (l-p)  yft  (l-p)Z  ft 

+  (l-p)  to  e  {(sine  +  y  cosc)[A( — ^ - )  -  A(— •)  j 


+  (y wsinc  -  cosc)[B( 


y  ftQ  (1— p ) ' 


ft 


o . 


-)  “  B(y±)  ] 


(4.24) 


where 


A(u)  = 


00  i/2  or  (i-l)/2  .  i-k  (i-2k)  i 

z  z  ^  L) - y - U 


i=0  k=0 


2k! (i  -  2k)!  2i+l 


and 


B  (u)  = 


i-1  i-2 

.  —  or  — 

Z  Z 

i=l  k=0 


,  - . i-k-1  i-2k-l  i 
(-1)  y  u 


(2k  +1)!  (i  -  2k  -  1)  !  (2i  +  1) 


In  the  summation  the  upper  limit  that  gives  an  integer  value  is 


used . 


' 


- 
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Given  the  nondimens ional  tension  x  then  normal  load  f  can 
be  found  from  equation  (4.12).  If  the  load  to  is  zero  then  the 
solution  is  much  simpler.  From  equation  (4.23)  the  tension  is  given 
by 


ft  -yft  p  (2-p)/2  ft 
/  O  \  o  o 

T  =  (t  -  — )  e  +  — 

o  2  p 


(4.25) 


and  from  equation  (3.12)  the  normal  load  is  then 


ft  -yft  p(2-p)/2  ft 

f  =  ft  (1-p )  [  (x  -  e  °  +  — ] 

o  o  2  y 


(4. 26) 


4.3.3.  Normal  Loading  (y  =  0) 


When  the  external  load  is  normal  to  the  rod  equation  (4.11) 
can  be  integrated  in  general  to  give 


x  =  C  +  —  (”)  +  o of  sin6  dp 

2  dp 


(4.27) 


for  the  nondimensional  tension.  The  integral  in  the  last  term  corresponds 
within  a  constant,  to  the  value  of  the  nondimensional  vertical  coordinate 
The  external  normal  force  f  is  given  by  equation  (4.12)  and  the  shear 
force  and  bending  moment  can  be  calculated  from  equation  (4.8). 
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4 . 4  The  General  Case. 

For  the  cases  in  which  the  shape  is  such  that  equations  (4.11) 
or  (4.13)  cannot  be  integrated  analytically  then  the  integration  must 
be  done  numerically.  Often  the  specified  shape  is  given  in  terms  of 
the  cartesian  coordinates  instead  of  the  function  0(s). 

When  the  shape  is  given  as  y  =  k(x)  a  Simpson ’s  rule  integration 
procedure  is  used  to  first  determine  the  length  of  the  rod  L  from  the 
equation 

d  1/2 

L  =  /  [1  +  (-j4  ]  dx  (4.28) 

o 

where  d  is  the  horizontal  distance  between  the  ends  of  the  required 
shape.  Once  the  length  is  known  the  shape  function  is  put  into  the 
form 


A 


(4.29) 


From  equation  (2.7)  in  terms  of  the  nondimensional  coordinates 


.dA 

0  =  arc  tan  (— ) 
d  4 


d£  0 

dp 


[1  +  ( 


dA 

d£ 


2  -1/2 

)  3 


and 


(4.30) 


< 
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Equations  (4. 30 )  can  be  used  to  obtain  the  following 
derivatives . 


de  =  (dX  2  -3/2 

dp  dc2  d? 


d20  d3X  .  ,dA  2  2  dX  ,d2X  dX  2 

[1  +  (^  !  -  3  «  (~2)[1  +  (d?>  ] 

d£ 


“3 


2  3 

dp  d£ 


d39  A  r  ,dA*2  ' 

,  3  ,c4  U  W  J 

dp  d£ 


'5/2 .  UAA  §  [i  +  A2] 

d£3  d£2  d£  dii 


-7/2 


-  3  (^-|)3[1  +  (%2]  '  +  lsA-|)3(-^)2[l  +  (~A2]  9/' 

d£2  dS  d?2  dC  d5 


(4.31) 


These  express ions can  be  substituted  in  equation  (4.13)  to 
give  an  expression  of  the  following  form  for  the  tension  t  as  a  function 

of  e 

y  0  ^ 

t(£)  =  e"y9(5)noe  °  +  /  R(£)  d£]  (4.32) 

o 

From  equation  (4.30)  the  coordinate  p  is  given 


£  Hi  9  1/2 

p(£)  =  /  [1  +  (-dj)  )  d£. 


(4.33) 


where  p  =  0  at  £  =  0. 
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The  integrals  in  equation  (4.32)  and  (4.33)  can  be  evaluated 

numerically  using  a  Simpson  rule  integration  procedure. 

The  normal  load  f(£)  can  then  be  evaluated  by  substituting 

equations  (4. 30)  and  (4.31)  into  equation  (4.12).  Similarly  the  moment 

ft(£)  and  shear  force  v (£)  can  be  evaluated  by  substituting  for  ■—  and 
9  dp 

(10 

— in  equation  (4.8)  from  equation  (4.31). 


dp 


4. 5  An  Example  Problem  with  Shape  Specified 

As  an  example  consider  the  family  of  shapes  shown  in  Figure  17. 
The  shape  is  given  by  a  seventh  order  polynomial 

7 


y  =  2  a .  x 

J  i 

i=0 


i 


(4.34) 


which  satisfies  the  conditions 


y (0)  =  y'(0)  =  y"  (0)  =  0 


y(d)  =  h  y'(d)  =  y"(d)  =  0 


(4.35) 


The  conditions  in  (4.35)  correspond  to  a  rod  that  has  zero 


g^Qpg  and  moment  on  the  two  ends.  Also  the  rod  has  a  change  in  elevation 
of  h  for  a  horizontal  distance  d.  The  second  last  condition  ensures 
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c  =  tan  50 
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FIGURE  17 :  Polynomial  Shapes 
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symmetry  about  the  mid  point  while  the  last  is  used  to  give  a  family 
of  shapes. 

Using  these  conditions  the  coefficients  in  (4.34)  are 


a  = 
o 


0 


.2 

a3d 

70 

-  32 

a3 
a,  d 

4 

-315 

160 

a4 

S5d 

= 

546 

-288 

a.d^ 

6 

-420 

224 

a-j&j 

120 

-  64 

h/d 

C 


(4.36) 


Given  Cs  h  and  d  the  length  of  the  rod  L  can  be  found  from 
equation  (4.28). 

The  nondimens ional  coordinates  of  the  shape  are  then  given 


by 


7 

A  =  Zb.  C1 
i=0  1 


where  b . 

i 


i  =  0,  1,  2,  ...  7. 


(4.37) 


Figure  17  shows  the  shapes  of  the  rod  in  nondimensional  form 


with  h/d  =0.4  for  three  values  of  C. 
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The  forces  required  to  maintain  these  shapes  can  be  calculated 

using  the  procedure  of  section  4.4.  For  this  problem  the  normal  load  is 

positive  and  negative  in  different  regions  of  the  rod.  To  keep  the 

frictional  force  applied  in  one  direction  along  the  rod  p  was  taken 

as  positive  in  the  region  where  the  normal  force  was  positive  and  as 

negative  where  the  normal  force  was  negative.  Figure  18  shows  the 

normal  forces  required  for  the  case  of  to  =  100  and  p  =  0.1.  The  initial 

tension  at  p  =  0  was  taken  as  x  =10.  The  tension  in  the  rod  for  three 

o 

values  of  C  is  shown  in  Figure  19.  The  shear  force  and  bending  moments 
along  the  rod  are  shown  respectively  in  Figures  20  and  21. 
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FIGURE  18:  Nondimensional  Normal  Loading  for  Polynomial  Shapes 
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FIGURE  19:  Nondimensional  Tension  for  Polynomial  Shapes 
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FIGURE  20:  Nondimensional  Shear  Force  for  Polynomial  Shapes 
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FIGURE  21:  Nondimensional  Moment  for  Polynomial  Shap 
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CHAPTER  V 

EXPERIMENTAL  CONSIDERATIONS 

5. 1  Experimental  Procedure 

Experimental  verification  of  the  numerical  solution  for 
the  stiffened  catenary  problem  considered  in  section  3.5.3.  is  given 
in  this  chapter. 

The  boundary  conditions  for  this  problem  can  be  duplicated 
experimentally.  The  procedure  followed  was  to  suspend  a  6  foot  long, 
1/16  inch  diameter  steel  rod  by  thin  threads  at  both  ends.  A  Cooke, 
Troughton  and  Simms  precision  level  was  used  to  position  the  ends  of 
the  rod  at  the  same  level.  The  coordinates  of  six  points  along  one 
half  the  rod  were  measured  using  a  cathetometer  (made  by  Gaertner 
Scientific  Corporation)  and  a  30  inch  long  brass  scale  clamped  to  the 
edge  of  a  2.5  foot  by  5.0  foot  surface  table  as  shown  in  Figure  22. 

The  cathetometer  was  used  to  measure  the  vertical  coordinates  of 
points  along  the  rod.  Its  telescope  and  micrometer  screw  have  a 
displacement  of  only  2  inches,  thus  larger  displacements  had  to  be 
measured  in  several  increments.  The  horizontal  coordinates  of  points 
along  the  rod  were  measured  using  the  brass  scale  in  conjunction  with 
a  vernier  scale  attached  to  the  base  of  the  cathetometer.  These 
coordinates  were  measured  for  three  different  slopes  of  the  supporting 
threads.  These  slopes  were  determined  by  measuring  the  horizontal  and 
vertical  distances  between  two  points  on  the  threads  using  the  brass 
scale  and  cathetometer.  The  experimental  data  is  given  in  Appendix  II. 
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FIGURE  22: 


Experimental  Apparatus 
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5 . 2  Calculations  and  Results 

The  nondimensional  coordinates  of  the  deformed  rod  are 

given  by 


X*  =  X/L 


Y*  =  Y/L 


(5.1) 


where  X,  Y  and  L  are  shown  in  Figure  11.  The  values  of  these 
coordinates  at  points  along  the  rod  as  calculated  from  the  experimental 
data  using  equation  (5.1)  are  shown  in  Figure  23. 

The  angle  a  defined  in  Figure  11  is  given  by 

a  =  arctan  (A  /A  )  (5.2) 

Y  A 


where  A  and  A  are  the  horizontal  and  vertical  distances  between  two 
X  Y 

points  on  the  supporting  threads.  The  values  of  a  calculated  from 
the  experimental  measurements  of  A  and  A  are  given  in  Appendix  III 
(Table  8). 

The  nondimensional  weight  of  the  rod  as  defined  in  section 
3.5.3.  is  given  by 


W* 


64  (WL)  L 


tt  D  E 


(5.3) 


for  a  circular  rod  of  diameter  D.  For  the  rod  used  in  the  experiment 


■ 
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with  E  =  29.8  x  10  psi,  the  nondimensional  weight  is  W*  =  1.805. 

Using  this  value  of  W*  and  the  values  of  a  given  in  Table  8 
the  numerical  solution  for  the  shapes  of  the  rod  can  be  evaluated. 

These  shapes  are  compared  to  the  experimentally  determined  points 
in  Figure  23. 

5. 3  Error  Analysis 

The  deformed  shape  of  the  rod  has  been  determined  experimentally 
by  measuring  the  coordinates  of  points  along  the  rod.  It  has  also  been 
found  by  experimentally  determining  the  loading  parameters  for  the  rod 
and  then  using  the  numerical  method  to  calculate  the  shape  of  the  rod. 

Both  solutions  are  based  on  experimentally  measured  quantities  which 
have  some  error  associated  with  them. 

If  x  is  the  measured  value  of  a  quantity  with  correct  value 
q,  the  error  is  defined  as 

e  =  |  q  -  x  |  (5.4) 

If  the  quantity  q  is  a  function  of  several  quantities,  say  q  =  F(a,b,c,d) 
then  the  error  e(q)  is  given  by 

o  d  2  2 

e(q)  =  Z  e  (n) 

n=a 


(5.5) 


\ 


- 


' 


when  the  errors  e(a),  e(b),  e(c)  and  e(d)  are  symmetric  and  independent 
(see  Reference  [5]). 


The  largest  differences  between  the  experimentally  determined 

coordinates  and  the  values  found  using  the  numerical  method  occur  at 

the  lowest  point  of  the  catenary.  The  values  of  the  X  coordinates  and 

vertical  deflection  at  this  point  are  h  and  d  as  shown  in  Figure  11. 

The  values  of  these  quantities  measured  experimentally  are  referred 

to  as  hr  and  d  ,  and  those  calculated  using  the  numerical  method  as 
b  b 


h„  and  d 


The  nondimens ional  values  of  these  quantities  are 


and 


(5.6) 


Using  equations  (5.5)  and  (5.6)  the  relative  errors  in  h*  and  d* 

b  b 


are  given  by 


£2  (h*)  =  s2(h)  +  s2(L) 


(5.7) 


£2(d*)  =  t2(d)  +  £2  (L) 

Ej 


where  the  relative  error  of  the  quantity  q  is  defined  as 


(5.8) 


v 


. 
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e(h),  e(d)  and  e(L)  are  the  relative  errors  in  the  measured  values  h, 
d  and  L. 

Determining  e(h*)  and  e(d*)  is  more  complicated.  The 

N  N 

quantities  h*  and  d*  are  functions  of  the  quantities  W*  and  a  which  in 
turn  are  functions  of  several  experimentally  measured  quantities. 

Using  equations  (5.3),  (5.5)  and  (5.8)  gives  the  relative 
error  in  the  nondimens ional  weight  as 

e2(W*)  =  e2 (WL)  +  4e2(L)  +  e2(E)  +  16e2(D)  (5.9) 


Similarly  the  relative  error  in  a  is  given  by 

r . 2  2 , .  .  J  2  2  .  .,1/2 
[A  e  (A  )  +  A  e  (A  )] 

e  (a)  =  ±  — - ^ - \ - - -  (5.10) 

(4 +  v  a 

After  calculating  e(a)  and  e(W*)  the  relative  errors  in  h*  and  d*  are 
given  by 


e  (d*)  = 


3d* 
a _ N 

d*  8  a 

L  N 


£  (a)  + 


d*  8W* 

L  N 


e2 (W*) 


c  (h*>  - 


8h* 
a _ N 

h*  8  a 
N 


e  (a)  + 


‘W*  _^N 

h*  8W* 
N 


e2(W*) 


(5.11) 


The  partial  derivative  must  be  evaluated  numerically  from 
the  numerical  solution.  The  values  of  h*,  d*  and  the  partial 
derivatives  calculated  with  the  three  values  of  a  arid  W*  =  1.805  are 


given  in  Appendix  III. 
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Using  the  data  in  Appendix  II  the  relative  errors  e(h*)  and 

E 

e(dg)  were  calculated  from  equation  (5.7).  These  are  given  in 
Table  5. 

The  relative  errors  e(W*)  and  e(a)  were  calculated  from 
equations  (5.9)  and  (5.10)  from  the  experimental  data  with  e(E)  =  ±  1%. 
The  relative  errors  e(h*)  and  e(d*)  were  then  calculated  from  equation 
(5.11)  using  the  quantities  in  Appendix  III.  These  relative  errors 
expressed  as  a  percentage  are  given  in  Table  5. 


TABLE  5 


Percentage  Errors  in  the  Numerical  and  Experimental  Solutions 


Case  1 

Case  2 

Case  3 

c(d*} 

±  0.5 

±  0.3 

±  0.2 

e(d*) 

±  0.6 

±  1.4 

±  1.6 

e<hf) 

±  0.1 

±  0.1 

±  0.1 

±0.01 

±  0.1 

±  0.2 

The  differences  between  the  numerical  and  experimental  values 
of  d*  and  h*  were  calculated.  In  all  cases  this  difference  was  less 
than  the  sum  of  errors  given  by  the  two  methods.  That  is 


dE  “  dN^  <  le(dE} I  +  'e(dN) 


V 


.. 


' 
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and 


|h*  -  h*|  <  |e(h*)|  +  |e(h*)| 

This  means  there  is  a  region  of  intersection  of  the  interval  d*  ±  e(d*) 

lL  ill 

and  the  interval  d*  ±  e(d*)  and  similarly  for  h*  and  h*.  It  can  be 

N  N  N  E 

concluded  that  within  the  accuracy  of  this  experiment  no  significant 
difference  was  found  between  the  numerical  solution  and  the  experimental 
results . 


* 
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CHAPTER  VI 

CONCLUDING  REMARKS 


6. 1  Summary 

The  pipeline  problem  mentioned  at  the  beginning  of  this 
thesis  requires  solution  to  two  types  of  problems  of  nonlinear  bend¬ 
ing  of  thin  rods.  These  have  been  considered  for  the  case  of  an 
initially  straight  rod. 

A  general  numerical  method  has  been  presented  to  solve  the 
problem  of  determining  the  deformed  shape  of  an  initially  straight 
rod  when  the  loading  is  given.  This  method  can  be  used  to  solve 
problems  of  rods  with  varying  flexural  stiffness  and  with  any 
combination  and  number  of  concentrated  forces,  couples  and  nonuniform 
distributed  loads  applied  to  the  rod.  Because  of  the  use  of  only 
one  segment  of  the  rod  at  a  time  in  the  numerical  method,  the 
calculations  can  be  done  on  a  small  computing  facility  such  as  the 
HP9820A  desk  top  programmable  calculator. 

The  problem  of  finding  the  distributed  loading  necessary 
to  deform  a  straight  rod  to  a  specified  shape  has  also  been  considered. 
Normal  and  tangential  distributed  loads  in  combination  with  a  uniform 
gravitational  load  were  used.  A  general  numerical  integration 
technique  was  developed  for  problems  where  the  specified  shape  is 
given  in  Cartesian  coordinates.  From  the  example  problem  considered 


t 


* 
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in  section  4.5  it  can  be  seen  that  small  changes  in  the  specified 
shape  of  the  rod  lead  to  large  changes  in  the  required  distributed 
loads . 

An  experimental  verification  of  the  numerical  method  for 
the  problem  of  the  stiffened  catenary  has  been  given.  The  experimentally 
determined  deformed  shape  and  the  numerical  solution  agree  within  the 
estimated  experimental  errors  for  the  particular  rod  considered. 

6. 2  Areas  of  Further  Research 

The  methods  and  solutions  mentioned  above  can  be  extended 
to  consider  rods  which  are  initially  curved.  For  the  problem  of 
finding  the  required  loading  to  bend  a  rod  to  a  specified  shape  it 
is  possible  to  consider  the  shape  given  in  polar  coordinates  or 
through  parametric  equations.  In  the  solution  of  this  problem  the 
function  defining  the  shape  of  the  rod  must  be  differentiated  several 
times.  Also  if  the  shape  is  specified  as  a  set  of  discrete  points 
then  some  interpolation  function  must  be  used  to  determine  the  required 
loading.  These  steps  can  be  avoided  by  using  a  variation  of  the 
numerical  method  used  to  solve  the  problem  of  unknown  deformed  shape. 

If  the  rod  is  divided  into  a  series  of  segments  each  loaded  by  a 
uniformly  distributed  or  concentrated  load  then  the  values  of  these 
loads  could  be  determined  using  a  one  dimensional  secant  method  so 
that  the  rod  passes  through  specified  points  of  the  required  shape 
function.  The  segment  solutions  needed  have  already  been  developed 
for  the  problem  of  unknown  deformed  shape. 


s 
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The  dynamic  problem  of  determining  the  natural  frequencies 
and  modes  of  vibration  of  the  rods  mentioned  above  could  also  be 
considered.  This  may  be  of  some  importance  in  the  laying  of  offshore 
pipelines  as  they  can  be  excited  by  ocean  currents. 
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APPENDIX  I 

GENERAL  PROGRAM  FOR  THE  NUMERICAL  METHOD  OF  CHAPTER  III 

This  program  has  been  set  up  to  solve  the  boundary  value 
problem  as  in  section  3.4.  It  can  be  used  to  solve  problems  with 
varying  distributed  loads  and  nonuniform  cross  sections  by  consider¬ 
ing  the  rod  as  a  series  of  segments  with  uniform  cross  sections  and 
uniform  distributed  loads  which  approximate  the  distributed  loading 
and  rod  cross  sections. 

The  first  part  of  the  program  solves  the  initial  value 
problem  as  described  in  section  3.4.1.  This  can  be  used  with  loads 
X  and  in  the  segment  or  with  tangential  loads  £  and  n*  The 
second  part  of  the  program  is  used  to  find  the  unknown  initial 
conditions  in  the  boundary  value  problem.  It  has  been  written  for 
the  case  of  one  and  two  unknown  initial  conditions.  The  case  for 
three  unknown  quantities  can  be  set  up  using  the  equations  in  section 
3.4.2. 

The  program  has  been  written  with  the  form  of  language 
used  in  the  9820A  Hewlett  Packard  programmable  desk  computer.  The 
register  numbers  (Rl,  R2,  ...  etc)  used  in  the  actual  program  have 
been  replaced  by  the  variable  names.  The  arrow  used  in  the  program 
statements  indicates  that  the  expression  to  the  left  of  the  arrow 
is  calculated  and  assigned  to  the  register  representing  the  quantity 
to  the  right  of  the  arrow.  In  the  first  part  of  the  program  the 


v 
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subscript  j  which  was  used  in  section  3.3  and  section  3.4  to  indicate 
.  th 

the  j  segment  has  been  dropped.  Line  5  and  6  are  used  to  assign 

known  and  unknown  initial  conditions  at  the  beginning  of  the  rod  for 

a  particular  boundary  value  problem.  Line  7  assigns  the  value  of  the 

ratio  of  segment  flexural  rigidity  to  the  reference  rigidity  as  defined 

in  equation  (3.18).  Line  10  assigns  values  to  the  constant  segment 

distributed  loads  Y-  and  ■  or  n .  and  Line  21  is  used  to  terminate 

J  J  J  J 

the  element  series  solution. 

For  one  unknown  initial  condition  the  one  dimensional  form 
of  equation  (3.36)  given  by  equation  (3.42)  must  be  solved  for  Z  which 
is  assigned  to  the  unknown  initial  condition  in  line  5  or  6.  The 
quantity  q*  is  the  specified  value  of  the  quantity  q  which  is  given 
at  the  end  of  the  rod.  One  of  the  quantities  B*,  D*,  G*,  X *,  Y*,  y 
at  the  end  of  the  rod  is  assigned  to  q  in  line  34  of  the  General 
Program  Part  2  (with  1-D  Secant).  The  one  dimensional  form  of  the 
iteration  given  by  equations  (3.37)  and  (3.38)  is 

-i+l  _  _i  _  (Z^Z1-1)  pCZ1) 

L  ~  L  i  i-1 

pCZ1)  -  P (z1  1) 

Two  initial  guesses  Z°  and  Z  must  be  used  to  start  this  iteration. 

With  two  unknown  initial  conditions  the  General  Program 
Part  2  (with  2-D  Secant  Method)  must  be  used  to  solve  the  equations 


p  (Z  ,  Z2)  =  q1(Z1,  Z2)  ~  q-*  =  0 


V 


#1 


*  • 


■ 
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P2(Zr  z2}  =  q2(Zl’  V  “  <*2*  =  0 


where  are  the  unknown  initial  conditions  at  the  start  of  the 

rod  assigned  in  lines  5  and  6  and  and  are  the  required  values 
of  q^  and  q^  which  are  the  quantities  specified  at  the  end  of  the  rod 
and  assigned  in  lines  34. 

The  iteration  used  in  this  program  is  the  two  dimensional 
form  of  that  given  by  equations  (3.37)  and  (3.38). 

z{i+1)  =  z/  -  [J22  ppz/,  z/)  -  J12  p^1.  z2i)]/D 

z2(i+1)  -  z/  -  [ju  p2(0  0  -  j21  px<0  0J/D 

where 


J 


J 


J 


J 


11 


12 


21 


22 


[ppz,1,  0 

Iff2/’  0 

[p2(z1i*  0 
tp2(z1i,  0 


-  P1(z1i  \  z2i)]/(z1i-z11  b 

-  Pp0  z2i-1)]/(z2i-z2i-1) 

-  P2(Z11"1,  z21)]/(z11-z11_1) 

-  P2(Z11,  z21_1)]/(z21-z21"1) 


and 


■ 


'.t 


■ 


'  -• 
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°  jll  j22  J21  j12 


Initial  guesses  (Z^  ,  )  and  (Z^  ,  Z^  )  must  be  used  to  start  this 

iteration. 


General  Program  -  Part  1  (with  distributed  segment  loads  Xj  and  ^ ) 

N  =  number  of  segments  used.  -  number  of  points  at  which 

data  is  printed. 


1. 

Enter.  N,  N  , 
P 

[loading  parameters  for 

a 

particular 

problem] . 

2. 

Print,  N,  N  , 

p 

[loading  parameters  for 

a 

particular 

problem] . 

3. 

GOTO  37. 

4. 

1  J. 

5. 

[  ]  ->  X*;  [ 

]->¥*;[  ]  -*»  y. 

6. 

[  ]  +  A*  [ 

]  ->  C*  [  ]  ->  F*. 

7. 

[  3  +  r. 

8. 

A .*/  (rN)  ->  a; 

C*/  (rN2)  ->  c;  F*/(rN2)  -> 

£ 

-L  . 

9. 

-  (3a  +  f ) /6 

3  • 

10. 

[  3  ->  x;  t 

3  V- 

11. 

4  ->  k. 

12. 

3  *  ak-4;  a  ^ 

ak-3;  f^2  ^  ak-2;  (ca 

X3 

~  /6  ->  ( 

Vl- 

13. 

ak~4  +  ak-3  + 

ak-2  +  Vl  +• 

14. 

ak-3  +  2ak-2 

+  3cx.  ,  ->■  b . 
k-1 

15. 

2\-2  +  6ak~l 

+  g. 

. 


■ 
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16'  ak-4  +  \-3/2  +  ak-2/3  +  ak-lM  •*  h' 

17 e  (cot^_2  “  X^k_ 3)/  (k(k-l))-*-  cik. 

18.  <p  +  -+  (f>;  b  +  ka^  ->  b. 

19.  g  +  k(lc-l)  ->  g;  h  +  a  /  (k+1)  ->  h. 

20'  “k-2  ^  \-3;  Vl  "  "k-2;  “k  "  Vl;  k+1  +  k' 

21.  If  ABS[ak  +  ak_1  +  «k_2]  >  [  ];  GOTO  17. 

22.  c  +  f 8  -  g(j)  ~  x  d. 

23.  X*  +  (cos(y-B)  -  hsin(y-3) )/N  ->  X*. 

24.  Y*  +  (sin(y-B)  +  hcos(y-3))/N  ->  Y*. 

25.  y  -  8  +  <f>  -*•  y. 

26.  brN  ->  B*;  drN2  ->  D*;  grN2  G*. 

27.  If  Q  i  1;  GOTO  31. 

28.  If  ABS  sin (  JttN^/N)  >  IE-6;  GOTO  31. 

29.  PRINT,  J,  B*,  D*,  G*. 

30.  PRINT,  X*,  Y*,  y. 

31.  B*  ->  A*;  D*  ->  C*;  G*  ->  F*. 

32.  If  J  i  N;  J  +  1  ->  J;  GOTO  7. 

33.  If  Q  =  1;  GOTO  [END]. 


General  Program  -  Part  1  (with  distributed  segment  normal  and  tangential 

loads  n .  and  £ . ) 

J  i 


Modifications  to  General  Program  -  Part  1  (with  loads and  ip  ) . 

10.  [  ]-*£;[  ]  n. 

12.  3  ->  ak_4?  a  +  ak_3;  f/2  *  \-2;  ( ca" n)/6  **■  \_i* 


• 

. 

. 


17. 


(c\-2  "  C  V3(1+1/k))/(k(k"1))  +  V 

22.  c  +  £3  -  g4>  -  £  -  nh  ->  d. 

General  Program  -  Part  2  (1-D  Secant  Method) 

34.  [  ]  ->  q. 

35.  If  R  =  1;  GOTO  41. 

36.  If  R  =  2;  GOTO  43. 

37.  ENTER  Z°,  Z1,  q*. 

38.  PRINT  Z°,  Z1,  q*. 

39.  Z°  ->  Z1”1;  Z1  Z1. 

40.  Z1_1  ->  Z;  1  ->  R;  GOTO  4. 

/i  .i-1 

41.  q  -  q*  ->  p 

42.  Z1  ->  Z;  2  R;  GOTO  4. 

•  «  • 

43.  q  -  q*  ->  p  ;  PRINT,  i,  p1,  Z1. 

44.  If  ABS  p1  <  [  ];  1  Q;  GOTO  4. 

45.  Z1  -  ((Z1-Z1"1)/(p1-p:L_1))  p1  +  Zi+1. 

,r  7i-l  i  i-1 

46.  Z  ->  Z  ;  p  p 

47.  Z1+1  ->  Z1. 

48.  1  +  i  ->  i;  GOTO  42. 

49.  END. 


General  Program  -  Part  2  with  Two-Dimensional  Secant  Method 


34.  [  ]  ->  q±;  [  ]  ■>  q2. 


35.  If  R=  1;  GOTO  44. 


* 


■ 


•  .  ■  .  ..  - 
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36.  If  R  =  2;  GOTO  48. 


37. 

38. 

39. 

40. 

41. 

42. 

43. 

44. 

45. 

46. 

47. 

48. 

49. 

50. 

51. 

52. 

53. 

54. 

55. 

56. 

57. 

58. 

59. 


If  R  =  3;  GOTO  50. 


ENTER,  Z  °,  Z  °,  Z  1 ,  Z  X,  q  *  q  * 

Z.  JL  4.  JL 

PRINT,  Z  °,  Z  °,  Z  1,  Z  1,  q  *  q  * 

j-  4.  -L  -L.  Zm 


1  i;  0  ->  Q. 

zi°  +  zi1_1;  z2°  z21"1- 

Z  ->  Z  1;  Z  ->  Z  1. 

d-  X  Z.  4. 

Zi1  ^  zi;  Z21  Z2;  1  R;  G0T0 

qi  “  qi*  **  p^1*1);  9 2  “  q2*  ^  P2^isi^ 

•  •  •  • 

PRINT,  i,  p  ^  p^,  Z11,  Z2X. 

•  c 

If  ABS[Pl  ]  <  [  ]  and  ABStp^]  <  [  ];  1  +  Q 

Z11  ->  Z  ;  Z21-1  ->  Z2;  2  ->  R;  GOTO  4. 
q  -  q  *  ->  p  (i,i-l);  q  -  q2*  ->  p  (i,i-l). 
Z11_1  ->  Zx;  Z21_1  -*  Z2;  3  ->  R;  GOTO  4. 

q^  -  q^*  p^(i-i,i);  q2  ~  q2*  p2(i-i,i). 

[Pl(i,i)  -  P1(i-I,i)]/(Z11-Z11'1)  +  JlV 
[p3  (i,i)  -  P1(i,i-1)]/(Z21-Z21"  )  J12. 


[p2(i,i)  -  P2(i-I,i)]/(Z11-Z11  x)  •>  J 


21' 


[p2(l,i)  -  P2(i)i-1)]/(Z21-Z21  1)  ->  J 


22' 


J11  J22  J21  J12  +  D‘ 

Z1  ~  [J22  Pl(i>i)  "  J12  P2(i’1)]/D  *  Zi 

Z2i  "  tJll  P2(i,i)  '  J21  Pi<i>i)]''D  -*  z2 
Z  1  -  Z  1_1;  Z  1  -  Z  1_1. 


i+1 


i+1 


;  GOTO  4. 


. 
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60.  i  +  1  +  i;  GOTO  43. 

61.  END. 


Specialized  Lines  to  Solve  Particular  Boundary  Value  Problems 

Some  lines  of  the  general  program  must  be  specialized  for 
solving  a  particular  problem.  These  lines  are  shown  below  for  the 
problems  treated  in  section  3.5. 


Elastica  Problem 


General  Program  - 

Part  1  (with  distributed 

1. 

ENTER  N,  N  ,  P*. 

2. 

PRINT  N,  N  ,  P*. 

5. 

0  -*  X*;  0  Y*;  Z 

y. 

6. 

0  ->  A*;  -  P*  cosy 

->  C*;  -  P*  sina  ->  F*. 

7. 

1  ->  r. 

10. 

0  ->  x»  0 

General  Program  -  Part  2  (with  1-D  secant  method). 


34. 


Y  ->  q. 


. 
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Cantilever  with  Uniform  Normal  Load 


General  Program  -  Part  1  (with  distributed  loads  £ ^  and  n ^ ) . 

1.  ENTER  N,  N  ,  P*. 

P 

2.  PRINT  W,  N  ,  P*. 

P 

5.  0  ->  X*;  0  -*  Y*;  Z  ->  y. 

6.  0  ->  A*;  0  ■>  C*;  0  F*. 

7.  1  ->  r. 

10.  0  +  K;  P*/(rN3)  ->  n. 

General  Program  -  Part  2  (with  1-D  secant  method). 

34.  y  ->  q. 

Stiffened  Catenary 

General  Program  -  Part  1  (with  distributed  loads  Xj  and  • 

1.  ENTER  N,  N  ,  a,  W*. 

2.  PRINT  N,  N  ,  a,  W*. 

5.  0  ->  X*;  0  -*•  Y* ;  Z  ->  y. 

6.  0  ->  A*;  W* (etna  cosy  -  siny)  +  C*;  W*(cosy  +  etna  siny)  -*■  F*. 

7.  1  ->  r. 

10.  -  W*  sin(y-3)/(rN3)  ->  x;  w*  cos  (y-3)/ (rN3)  -> 


107. 


General  Program  -  Part  2  (with  1-D  Secant). 
34.  y  -*  q. 


Nonuniform  Cantilever  (with  specified  distributed  load  and  concentrated 

end  load) . 


General  Program  -  Part  1  (with  distributed  loads  Xj  and  40 


1. 

ENTER 

N, 

N.  , 
P 

V* 

P*, 

a,  t . 

2. 

PRINT 

N, 

N  , 
P 

V* 

P*, 

a ,  t . 

5. 

0 

-> 

X* 

• 

y 

0 

Y*;  Z 

->  y 

• 

6. 

0 

-> 

A* 

1  • 

P*  cosot 

C*; 

_  p* 

7. 

t 

+ 

(J 

**“  • 

5)  (1 

-t)/N 

->  r 

• 

10. 

— 

W  * 

i-j 

sin  (y~ 

6)/N3 

+  X 

•  W  * 

’  WL 

-> 


4> 


General  Program  -  Part  2  (with  1-D  Secant) 


34.  y  ->  q. 


Nonuniform  Cantilever  (with  Yq  specified  and  no  concentrated  end  load) . 

General  Program  -  Part  1  (with  distributed  loads  \  and  4^)* 

1.  ENTER  N,  N  ,  YQ »  t. 

2.  PRINT  W,  N  ,  Y  s  t. 

P  o 

5.  0  ->  X*;  0  ->  Y*;  YQ  Y* 

0  ->  A*;  0  ->  C*;  0  ->  F*. 


6. 


■■ 
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7.  t  +  (J-.5)  (l-t)/N  ->  r. 

10.  -  Z  sin(y-3)/N^  ->  x  1  Z  cos(y~3)/N^  ->  ip. 


General  Program  -  Part  2  (with  1-D  secant) . 


34.  y  ->  q. 


•  - 


109. 


APPENDIX  II 

EXPERIMENTAL  DATA 

The  rod  used  in  the  experiment  is  composed  of  two  1/16  inch 
diameter  welding  rods  welded  together  at  the  mid  point  of  the  rod. 
Data  for  this  rod  is  given  below: 

weight  of  half  rod  WL  =  28.2  ±  0.1  grams 
length  of  half  rod  L  =  36.0  ±  0.03  inches 
diameter  of  rod  D  =  0.0625  ±  0.0002  inches. 

To  determine  the  slope  of  the  threads  supporting  the  rod 
the  horizontal  and  vertical  distance  A  and  A  between  two  points  on 

A  I 

the  thread  were  measured.  These  distances  are  given  in  Table  6. 

The  coordinates  of  the  deformed  shape  of  the  rod  measured 
for  the  slopes  of  the  threads  given  by  Table  6  are  shown  in  Table  7. 


TABLE  6 

Horizontal  and  Vertical  Distances  Between  Points  on  Thread 


Case 

A  (±  0.02  in. ) 

A 

A.y(±  0.002  in.) 

1 

3. 66 

1.040 

2 

1.55 

1.900 

3 

0.00 

1.241 

. 

. 


- 


' 


. 

Measured  Coordinates  of  the  Deflected  Stiffened  Catenary 
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APPENDIX  III 

NUMERICAL  SOLUTION 

The  values  given  in  Table  8  were  calculated  using  the 
numerical  method  to  solve  the  stiffened  catenary  problem  of  section 
3.5.3.  This  solution  is  compared  to  an  experimental  solution  in 
Chapter  V. 


TABLE  8 

Quantities  Calculated  using  the  Numerical 
Solution  for  the  Stiffened  Catenary  Problem 

with  W*  =  1.805 


Case  1 

Case  2 

Case  3 

a 

* 

0.2769 

0.8805 

1.571 

d* 

0.1033 

0.2219 

0.3259 

3d* 

3W* 

0.015 

0.069 

0.14 

3d* 

0.28 

3a 

0.15 

0.18 

h* 

0.9932 

0.9687 

0.9313 

3h* 

3W* 

-0.0020 

-0.020 

-0.060 

3h* 

3a 

-0.035 

-0.043 

-0.076 

* 


